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1. Introduction 


Economic problems characterizing nuclear power stations constitute the 
most heavily discussed themes in connections with the utilization of the new 
energy source. The discussion is still mainly of a theoretical character, no 
practical operational experience being available, which would permit the 
elaboration of rational economic parameters and the establishment of their 
interrelations for different types of nuclear reactors and nuclear power stations, 
as well as for conventional power plants. Uncertainties extend to both the 
domains of capital and to running costs. Among the problems discussed the 
valuation of plutonium produced in power reactors and its economic equiva- 
lent, i.e. plutonium credit stands at first place, at least so far as the uncer- 
tainty of valuation is concerned. In the present paper it is proposed to discuss 
this question, being of basic significance for the operation economies of nuclear 
power stations, but not yet cleared to its proper depth, only on theoretical 
foundations, according to the present state of affairs. 

The relevant literature [1] reflects different viewpoints, so far as the 
valuation of the plutonium produced in nuclear power stations is being con- 
cerned and these may be essentially grouped into the following classes : 

a) Most frequently the value of plutonium is referred to the fixed price 
of the isotope U2, thus taking a unit price for the plutonium at $ 15 000 to 
30 000 per kg. This unit price, however, — beyond ignoring the differences 
between the two fissile materials, so far as their equivalent heat content and 
their behaviour in reactors is concerned — should obviously be regarded as 
the upper limit of plutonium valuation, e.g. because it takes the plutonium 
yield — which in its pure state may become usable only after years in breeder 
reactors — at once at its full value, at the very moment of its production. 
A similar effect can be noticed as a consequence of this fact that the price 
of the isotope U?® is to-day determined by other considerations not of 
economic character. 

b) In contrast to the preceding conception, another group tends to 
find a relation between the price of plutonium and that of natural uranium 
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(of its fissile material content), in which case the price of natural uranium at 
$ 40 per kg will entail the value of plutonium per kg to be about $ 5600. 
This method of evaluation — beyond all the backdraws of the former method 
is also valid for the present case — does not take into account that the pluto- 
nium thus produced is a highly concentrated fissile material. In this case 
plutonium will obviously be underpriced. 

c) If one has a free hand to choose the appropriate energy carrier and 
the problem is being considered on a purely energetic basis, the price of pluto- 
nium may be referred as that of the cheapest fuel to be fired in a power station 
(e.g. fuel oil). If one wishes to take into consideration the fact that capital 
costs of present-day nuclear power stations usually amount to two or three 
times the specific investment necessary for conventional power stations, and 
therefore fissile material is being taken with one half of the usual value as a 
cost factor, then along this line — i.e. using only the concept of thermal 
equivalence — a plutonium price of $ 14 000 per kg wil) be obtained. 

d) Essentially a modification of the preceding paragraph can be seen 
in the method of evaluation, where the price of plutonium is determined on 
the basis of the lowest unit price of electrical power of a given country or a 
region. This procedure leads to the evaluating method proposed by the authors, 
if the quantity of plutonium to be evaluated is converted into electric power, 
not only by taking into account its theoretical best content, but also by con- 
sidering the reactor cycle utilized in the process. 

Evaluating procedures hitherto described are hampered by the follow- 
ing principal errors : 

— the evaluating process is based on more or less arbitrarily assumed 

energy carrier prices ; 

— the procedures ignore the fact that the behaviour of plutonium in a 

reactor differs from that of uranium ; 

— the isotope content of the plutonium produced, depends on the design 

and operation method of the reactor producing it ; 

— plutonium production is attached to significant capital and running 

costs ; 

— the plutonium produced — depending on the cycle to be chosen — 

many times will be utilized only after a substantial elapsed time etc. 

If, however, one defines as a leading principle that the plutonium produc- 
ed in nuclear power reactors is worth the value of the electrical power to be 
developed from it, then obviously the possibility of peaceful utilization will 
form the basis of evaluation and — as to be seen from our discussion — the 
grave uncertainty factors inherent with the preceding valuation procedures 
will greatly be eliminated. While difficulties are thus transferred to the fields 
of reactor engineering calculations and operational uncertainty factors, these, 
however, may be successfully approximated by mathematical and physical 
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methods and, therefore, resulting uncertainties will be less than with most of 
_ the evaluating procedures known from literature. For a similar evaluating 
method the value of plutonium is defined essentially by the way plutonium 
is utilized. Thus for determining the value of plutonium the method of utili- 
zation (i. e. the fuel cycle) must first be defined. With due regard to the large 
' number of theoretically available variations, it is proposed to base our investi- 
gations on the following assumptions — not necessarily the most favourable 
ones, so far as plutonium valuation is being concerned : 

— a heterogeneous thermal reactor is used for plutonium production ; 

— plutonium shall be reutilized in a similar reactor ; 

— in our reactor engineering calculations the simplifications to be 
enumerated in Section 3 will be introduced ; 

— in our computations only the varations of fuel costs will be taken 

- into account, while neglecting the effect of capital cost, because 
of lack of reliable data. (According to the authors’ opinion consider- 
ation of the latter cannot decisively influence the results to be obtained 
by these procedures.) 

— the fact that burnup levels in a reactor are frequently determined, 
not by reactivity requirements, but by corrosion considerations, will 
be ignored for the time being (the latter especially applics to fast 
reactors). 

Having introduced the above simplifying assumptions, two procedures 
will be analyzed later on, in order to present the principles of plutonium 
valuation : 

1. By means of reactivity analysis it is proposed to determine the burnup 
or irradiation level to be attained, in a reactor using natural uranium, and for 
given initial conditions. According to our assumption to be later motivated, 
the surplus burnup to be realized subsequent to the first burnup and the extrac- 
tion of fission products can be taken as net profit due to the plutonium produced 
in the reactor (para. 4.1.). 

2. Plutonium produced in a reactor burning natural uranium will be 
utilized — after extraction and mixing with natural uranium — as the basic 
fuel for a second fuel cycle. Thus a self-sustaining cycle may be realized 


(para. 4.2.). 
2. Method of plutonium valuation based on the fuel costs of electric power 
to be developed in a nuclear power station 
The specific variable cost resulting from the fuel consumption needed 
for producing the electric power fed by the nuzlear power station into the 


grid can be expressed by the following general relation : 
4170 5 : 
k= ae [pi + ty Ps + t2°Pr— Spa Pru + Pa] mills/kWha (1) 
27 


1* 
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In the above formula the following symbols have been adopted : 

Q — the burnup or irradiation level, i.e. the amount of heat energy to be 
produced from one ton of fuel in the reactor (MWd/ton) ; 

n — the efficiency of the nuclear power station (ratio of heat energy 
introduced into the heat exchanger to electric power fed into the 
mains, per cent) ; 

p, — fuel unit price taking reactor purity grade metal ($ per kg) ; 

i, — the number of fuel element production stages for realizing the burnup 
level Q; 

Pz — fuel element production costs (8 per kg) ; 

i, — number of reprocessing operations necessary for realizing a burnup 

level Q; 

Pr — costs of a single reprocessing of one kg fuel (8 per kg) ; 

gp, — plutonium yield of one kg nuclear fuel after reaching the burnup 
level Q (kg per kg) ; 

Ppu — value of 1 kg of plutonium produced in the burnt-up fuel ($ per kg) ; 

Pa — costs of extracting and disposing the radioactive wastes referring 
to one kg of fuel ($ per kg, the latter factor being, however, neglected 
during our investigations). 

Numerical examples as published in literature usually take a relation 
similar to that given above as a base for plutonium valuation, but the price 
of plutonium (pp,) is determined according to different assumptions as is 
enumerated under a) to d) in Section 1. In accordance with our aims, we are 
trying to find — in place of these — a different valuating process, where 
the value gp, of plutonium produced is defined by the electric power or excess 
heat energy (Q, MWd per ton) to be developed from it. It can be shown that 
plutonium value is thus defined by the relation 


Lei A 
Ppu = aa bs (py + Ps) — (41, “Ps Aig: pp) $/kg (2) 
pu 2 


where the new symbols denote : 


Ai, the number of additional fuel processing stages in order to realize 

the excess burnup level AQ to be obtained by utilizing plutonium 

Ai, the number of reprocessing operations needed for same. 

In our equation (2) giving the value of the plutonium produced we find, 
besides the terms Ai, and Ai, as defined above, the original burnup level of 
the reactor (Q), the excess burnup level to be realized from utilizing pluto- 
nium (AQ), the plutonium yield (gp,), as well as the costs of fuel (p;), fuel ele- 
ment production (p,) and reprocessing (pr). In the course of economic analysis, 
the following numerical values have been adopted for different cost factors, 
in accordance with literature sources [2.03] 2 
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Pi = *40 per kg (natural uranium) 

P, = %16 per kg (natural uranium) 

P, = %22 per kg (plutonium-enriched fuel) 

Pr = $13—25—37 per kg. 

In our proposed method of valuation the price of plutonium will be 
rationally determined as a function of reprocessing costs, the latter varying 
between wide limits. 

The other factors (Q, 4Q, gpy) are greatly dependent on the actual type 
of nuclear reactor, its operating cycle and the method of utilization of the 
plutonium produced — taking the simplifying assumptions already enumerated 
- into account. It seems therefore advisable to discuss — before turning to 
actual examples of plutonium valuation (see Section 4) — a few important 
reactor engineering problems necessary for the proper valuation of plutonium, 
namely: general analyses of plutonium production and determination of 
attainable burnup level (Section 3). 


3. Analysis of plutonium production in a nuclear reactor using natural 
or slightly enriched uranium 


3.1. Variations of fuel composition during the burnup cycle 


During permanent irradiation of reactor fuel elements several important 
changes occur in the fuel and in the interest of taking these quantitatively 
into account, the following simplifying assumptions must be introduced : 

a) Only thermal neutron flux (®) will be considered. This simplifying 
assumption, however, does not apply to several terms of the multiplication 
constant (k), thus the fast fission factor (¢), the resonance escape probability 
(p) and the neutron leakage factor (P). 

b) Variations of thermal neutron flux (®) within the reactor itself will 
be ignored, utilizing only the mean value. It is also assumed that the reactor 
is operated at constant heat release rate which obviously involves a neutron 
flux only slightly varying with time. Time dependence of the neutron flux 
has been, however, considered, as far as possible, when determining the burnup 
level. 

c) The quantity of the isotope U** placed into the reactor has a constant 
value and is identical to the starting value (i.e. s = const). 

d) Pu?8® is produced from U only as a result of the fission of U”* 
and Pu®®, 

e) Among the factors figuring in the multiplication constant €, p and 
P during irradiation remain unchanged. 

f) Fission products are produced according to the simplifying assump- 
tions as enumerated in para. 3.2. and their behaviour during irradiation is in 
conformation to the description given there. 
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It is known that in reactors burning natural or slightly enriched uranium 
as a fuel, plutonium production is a result of the so-called internal conversion 
and may be computed from several phases. Production of the isotope Pu?®9 
(its quantity being z) proceeds along the following lines : 

1. Thermal neutrons are captured by the isotope U** (the probability 
of capture being given by o,) and the isotope Pu? will be produced as a result 
of radioactive decay. The quantity to be produced within a differential time 
interval (dt) will thus be given by the term 


$+Ogg:P-dt 


where s = const, denoting the number of U®* atoms. 

2. During slowing down of the fission neutrons one part (1 — p) of 
them suffers resonance capture in the U8 according to the definition of the 
probability (p) of resonance escape. As fission neutrons represent products 
of the actual U2%5 (u) and the Pu? (z) content and of fast fission products (¢) 
(ignoring the fission of other materials, e. g. Pu™!), the number of Pu®® atoms 
produced by resonance capture. during the time element may be given by the 
formula : 


[> Og5°Ns+ 2*Fq9°%gl-P-e-[1 — p]-P-dt 


where 7 denotes the number of fission neutrons per captured neutron and P 
given the leakage factor. 

3. The actual number of Pu? atoms will be burnt up by the capture 
of thermal neutrons (0,9), partly through fission, partly through conversion 
into isotope Pu. This term is given by 


= Bit Cig Pr dt 
The approximate differential equation of Pu®® production will be therefore 
given by 
dz=s-0,g:P-dt + (u-045-3+2+Og9°N)-e-(1 — p)-P-D-dt — z- Oa9°P-dt (3) 


In view of the fact that the actual quantity of the U* isotope is defined by 
the relation 

WS Ug esas Poe (4) 
ve solution of the general differential equation may be given by the following 
relation : 


525° e&(1—p) Pry og, Pt — og,Pt ae 


&- (1 — p) 1P ang: Uy * Ogg [eas Pt @° (lap) Pity Gan Par con] 


- 


Gog — &+(1 — p)+ P+ 19g Gag — Soe + (5) 


t a as =X, [1 ef —p) P ny cay ® t—og, a | 


Gag —E> (lp): Pease sy 
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where the following notations have been introduced : 
U,y — the number of U* atoms present in the reactor at the start-up 
(equals unity, if not otherwise specified) ; 
% — the number of Pu*®® atoms present in the reactor at the start-up 
(equals zero, if not otherwise specified). 
In economic and operational research concerning reactor economy the so- 
called heavier isotopes of plutonium may play an important role, consequently 
their quantitative variations must also be taken into account. Differential 
equations established on the basis of nuclear physical relations as well as 
analytical or numerical solutions for these will be omitted here, only the final 
results are presented on Fig. 1, showing variations of Pu, Pu, Pu! and 
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Fig. 1. Variation of specific plutonium isotope yields versus irradiation level (for natural 
uranium ; reactor characteristics according to Appendix) 


Pu? isotope yields — referred to one initial U23> atom — as a function of 
burnup (®.t n per kilobarn). (Nuclear constants used for numerical solutions 
of the equations, as well as characteristics of the hypothetical reactor may be 
found in para. a, and b, of the Appendix.) 


3.2. Determinations of attainable burnup level 


The burnup or irradiation level to be realized in a reactor (in MWd 
per ton) is of decisive importance, so far as the unit cost of electric power 
produced by a nuclear power station and the utilization of fuel are concerned 
{see formula (1)]. Its actual value is defined — apart from variations in fissile 
materials as defined in the preceding Section — by time variations of the 
neutron — absorbing cross sections, and obviously by running requirements 
concerning reactor operation. As already pointed out, burnup level limitations 
due to corrosion are to be ignored. 

Laws governing the production — during a fission reaction — of fission 
products and of non-fissile isotopes of fissile elements will not be dealt with 
here. The problem is amply discussed by literature[4]. It should be, however, 
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noted here that reactor poisons — essentially defining burnup level — are 
usually classified into three group [5] 5 . 

a) The first group comprises all strong poisons, with an rove neutron 
capture cross section exceeding 1000 barns which, however, will become 
saturated within a short interval as a consequence of their nee neutron 
capture activity. Within this group, as usual, let us stress the importance of 
Xe! (gaxe 4 3,5-108 barn) and Sm“? (Gqsm & 5,3° 10° barn). 

b) The second group contains poisons of medium intensity (rise ogee 
10 barn), characterized by their way of formation and the compositions of the 


initial material. They cause, as a rule, a maximum poisoning effect a short 
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Fig. 2. Neutron capture cross sections of reactor poisons as a function of irradiation (natural 
uranium) 


while after starting reactor operation and this may decrease, somewhat later. 
Certain poisons, however, show a saturation effect. 

c) The third group comprises weak poisons (og < 10 barn), which do 
not show any saturation effects and which have a poisoning effect increasing 
with time. 

In our calculations the poisoning of the reactor has been determined 
on the basis of literature data [6]. Fig. 2 shows the variation of neutron 
capture cross sections characteristic of reactor poison groups and referring 
to an initial U5 atom. 

It is known that a nuclear reactor remains serviceable until its effective 


multiplication constant (kj = «-p-f:n- P) exceeds unity, that is until 
it has an effective excess reactivity, i. e. until 


ox = kery == ] > 0 (6) 


ANALYTICAL AND ECONOMIC INVESTIGATION OF PLUTONIUM PRODUCTION 271 


The overall number of fission neutrons for the reactor as a whole is 
given by the product 


prod = > (D1 -v) (7) 


where ©; denotes the macroscopic fission cross section of all fissile fuel ele- 
ments and y denotes the number of neutrons produced by a thermal fission. 
The sum of these products is called the thermal productive cross section of the 
reactor [7]. 
Let us denote the total absorbing cross section of fuel by © L750 that 
final factor 7 (i.e. the number of fission neutrons for a neutron captured by 
the fuel) for the reactor as a whole may he given by 


>Y 
eae (8) 
— a a,f 


According to its definitions, the thermal utilization factor (f) is given 
by the ratio of thermal neutrons absorbed by the fuel to the total number of 
neutrons absorbed by the reactor. Apart from the fuel (denoted by index f) 
neutrons are also being absorbed by reactor poisons (index t), by the moderator 
(index m), as well as by other structural elements, such as cooling agent, fuel 
fans, control rods, etc. (index r). Without considering the flux differences 
always present in the fuel elements and the moderator as well as other mate- 
rials, let us write down : 


f=s 
on! 


P Dass (9) 
7 7 
(D'a,F =e Da, i aL Da, m oy aa,r) 
The expression in the denominator of the preceding relation represents the 
overall absorbing cross section of the reactor as (Zaps) a whole, so that we 
have 


f= a Daf (10) 


as 
a abs 


Using the symbols as deduced above, the formula of the effective multiplica- 
tion constant may be written in another form : 


ae Bx fred , 
heyy =o: p-f-9-P =e: p-P AE (11) 
«abs 
and so for the effective excess reactivity we have : 
g, OPP: Arma = Lars (12) 


ays 
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The expression in the numerator represents the actual excess cross section 
of the reactor variable with time (2,,-) and therefore excess reactivity will 
be given by: 
y 
dere (13) 


ps abs 


The reactor will remain workable, until it still possesses an excess cross section, 
i.e. until its productive cross section exceeds the absorbing cross section. 
Particulars of the productive cross sections can be calculated by using the 
analysis given in Section 3.1. — at least with a fair degree of approximation — 
while calculations of the absorbing cross sections of an actual nuclear reactor 
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Fig. 3. Variation of producti i n n n 
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may also be undertaken if one takes into account the rate of poisoning — 
variable with time —, the variations of temperatures, steam generation within 
the reactor, etc., i.e. all operational factors, so that for a given case variations 
of the productive and absorbing cross sections, as well as of excess reactivity 
can easily be plotted as functions of time, resp. irradiation, and thus the burnu 
level to be realized within the actual reactor can be determined. ‘ 
A practical application of the preceding theoretical discussions is shown 


by Fig. 3, for a reactor complying with particulars of para. b, of the Appendix, 
where absorbing cross sections for poisons and non-poisons have been separa- 
tely defined. It can be read from this diagram that an irradiation level of 


a es 
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0,84 n per kilobarn may be realized with the specified reactor. Fig. 4 shows 
the variations of excess cross-section (0,,;) and excess reactivity (0,) with 
reference to one initial U?* atom and using the same data as in Fig. 3. 

As the procedure given in Section 3.1. permits the determination of 
changes in fuel composition for any arbitrary irradiation level, the appropriate 
integral curves will deliver the isotope quantities having suffered fission. By 
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Fig. 4, Variations of excess cross section and excess reactivity versus irradiation level 
(for natural uranium, reactor characteristics according to para. b, Appendix) 


knowing this, as well as the energy released during the fission of one gram of 
the given isotope (see para. c, Appendix), the heat energy produced can be 
determined. For our actual example an irradiation of 0,84 n/kb will yield 


some Q = 3350 MWd heat energy for one ton of natural uranium, according 
to the following distribution : 


fission of U**° 2612 MWd/ton 13°95 
fission of Pu? 708 MWd/ton PANES bs 
fission of Pu*4 30 MWd/ton 0,9% 


4. Plutonium valuation for different fuel cycles 


The governing principle of plutonium valuation as based on the concept 


of equivalent electric power has been described in Section 1, while calculation 
procedures were dealt with in Section 2. In the following paragraphs the valuat- 
ing procedure will be applied to two actual plutonium utilizing processes 


mentioned in literature, using the result of reactor computations given in 
Section 3. 
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4.1. Plutonium valuation of a fuel cycle with prolonged burnup 


After having been refined from accumulated fission products, reactor 
fuel becomes again serviceable. Using the particulars as well as the results 
of our preceding investigation of actual reactor service life it is found that 
subsequent to a fuel burnup of Q = 3350 MW4d/ton and after reprocessing the 
fuel from reactor poisons (0; * 25 barn) and refilling it to the reactor it will 
show an excess reactivity of 


8, = = = 0,019 


Fabs 
in the interest of a prolonged burnup period. 


Assuming that the fuel purified from poisons is replaced to the same 
reactor for another burnup period, as well as that the plutonium produced 
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Fig. 5. Variations of productive and absorbing cross sections versus irradiation level for 
a prolonged fuel cycle (for natural uranium, reactor characteristics according to para. b, 
Appendix) 


will not be extracted from the fuel, variations of the productive and absorbing 
cross sections for the additional burnup period can be determined in a way 
analogue to the analysis given in Section 3. It can be seen from Fig. 5 that in 
our example this used, but poison-free fuel, permits an additional irradiation 
of some 0,3 n/kb. Thus, one ton of natural uranium will yield an additional 
heat energy of 4Q = 1100 MWd. 

Until the extraction of poisons, i. e. up to an irradiation level of 0,84 n/kb 
or a burnup level of Q@ = 3350 MWd/ton, gp, = 2,4 - 1073 kg/kg plutonium 
will have been produced in the reactor. The additional heat energy of AQ = 
= 1100 MWd/ton as computed before, is essentially due to this quantity 
of plutonium. We may namely — with about the same reprocessing costs as 
necessary for extracting the poisons from reactor fuel, extract plutonium itself 
and sell or utilize it in another cycle (e. g. according to para. 4.2.). If, however 
plutonium is retained in the reactor fuel while extracting other fission prodncu 
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thus enabling the reactor to be operated so as to achieve an additional burnup, 
the former considerations will retain their validity. 

Through applying relation (2), plutonium value in the burnt-up fuel 
will be given according to the following table, for different reprocessing 
costs — taking into account the higher production costs of fuel elements 
mixed with plutonium. 


Preprocessing cost (pr) $/kg 13 25 37 
Plutonium value (pPu) $/kg —7530 —12340 —17350 


Results are shown by curve a on Fig. 6. As seen from it, plutonium utili- 
zation by means of a prolonged fuel cycle — for the present reactor — proves 
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Fig. 6. Plutonium value for different valuating processes as a function of reprocessing costs 


to be uneconomical even in case of low reprocessing costs, mainly because of 
inereased fuel element production expenses. It appears, therefore, to be logi- 
cal why this method of plutonium utilization was received unfavourably by 
the literature — at least for the time being [8]. 


4.2. Plutonium valuation for an equilibrium state fuel cycle 


A possible way of utilizing plutonium for energy production is the reali- 
zation of an equilibrium fuel system [9]. The governing principle in this system 
may be explained as follows: the reactor is started up with a fuel having 
a composition which ensures that — after reaching the burnup level — the 
plutonium produced in the reactor and mixed with natural uranium provides 
a fuel equivalent to the original fuel, so far as reactor operations are concerned. 
Thus the burnup of natural uranium in the reactor may be increased and also 
several subsequent, nearly equal running periods may be achieved. 
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With quite a degree of simplification, two reactor fuels (A and B) may 
be regarded as nearly equivalent, so far as reactor operation is being concerned, 
if — filled into a similar reactor — the effective multiplication constant of 
the reactor remains the same for both cases. (This identity from the point 
of view of reactor engineering does not necessarily entail that the two fuel 
types are equivalent so far as energy production is concerned.) Thus we have : 


keppa = Ea *PA‘Na -Pa-fa = hepa =€3* PBB -Pz-fe (14) 


If it is wanted to determine how many kilogramms of plutonium must 
be added to one ton of natural uranium in order to produce a fuel equivalent 
to slightly enriched uranium, then it may be assumed that 


Ex, =&33 Pa=Ppg; Pa=Pp 
Thus for the equivalence of reactor fuels we have the following condition : 


fans =fe-Ne (15) 


Using our previous notations, we may write down 


Ua'F,5 1+ Sa° Sag : N5°UA*Fa5 (16) 


fa na = 
Ua°Sgs + Sa°Sagt Fo ua: Oggt Sa *Sqp 


or 


Up * Fq5 + SB * Fag + 22+ Oa9 +28: Fao + ZB ° Oo + zh * Ga 


fa-Ns . 
WSs + $B*%qg + 2B Fag + 2R°Fq9 + 25°F + BBG gg + F 


. . Syste y 
5° Up Fs + Ny* 28° Fag + My 2B * Far 


UB Ogg + $B° Ogg + 2+ Faq + ZB Faq + 2B °F q1 + 2B ° Fa9 


(17) 
where 


F denotes the total capture cross section of neutron-absorbing materials 
in the reactor — except fuel — (assumed identical for both fuel types) ; 
index A: enriched uranium 
index B: mixed fuel consisting of natural uranium and plutonium. 

Results of the numerical solution of the above equation are shown in Fig. 7 

(based on particulars given in the Appendix). It may be inferred from this 
that the addition of about 5 kgs of Pu®®® isotope to one ton of natural uranium 
will result in a fuel which is equivalent to uranium enriched to 1,4 per cent, 
so far as reactor operational considerations are taken into account. (In order 
to simplify our calculations, the fact of plutonium consisting of a mixture of 
different isotopes has been ignored and only the isotope Pu8® has been taken 


into account when computing the equivalence of fuels. The error thus committed 
is, however, less than 1 per cent.) 
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The proper task of equivalence investigations is to establish the running 
time necessary for achieving an equilibrium fuel system, i. e. to determine how 
long a reactor must be kept in operation in order to obtain a quantity of 
plutonium which 

a) mixed with natural uranium provides an equivalent fuel for reactors 
started-up with enriched uranium fuel ; 

b) in case of reactors started-up with a mixed fuel of natural uranium 
and plutonium equals the quantity of plutonium filled-in prior at start-up. 
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Fig. 7. Plutonium content of natural uranium equivalent to slightly enriched uranium 


(For similar investigations losses occurring during plutonium extraction as 
well as fluctuations in plutonium composition may be neglected.) 

Solutions of the differential equations determining the quantities of 
isotopes — provided that the reactor already contains plutonium (z)) during 
start-up — have been given in para. 3.1. Let z denote the quantity of Pu®? 
in formula (5), assuming that there has been no plutonium present during 
start-up. In this case equation (5) may be transformed into the following 


general form : 
; a e&(1—Pp)P ty: Fay" P+t—oapPt | x (18) 


In a reactor using equilibrium fuel, we must have necessarily that subsequent 
to an irradiation @-t 

a Zoe 
This permits the computation of the quantity of plutonium necessary for the 
start-up : 


ik re - (19) 


ihe ef(1—P)PnctasPt—Fas «@.-t 


In the numerical example elaborated to illustrate the present method 
it has again been assumed that in the reactor — as specified by para. b), of 
the Appendix — the plutonium gp, = 2,4: 10~° kg Pu per kg natural U) 
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produced after the irradiation (Q = 3350 MWd per ton) determined during 
the service life investigations (Section 3.2) is mixed with natural uranium and 
replaced to the same reactor for repeated burnup. Variations of the productive 
and absorbing cross-sections of the reactor are shown in Fig. 8 for the reactor 
started-up with a mixture of plutonium and natural uranium. Maximum 
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Fig. 8. Variation of the cross-sections of an equilibrium system of plutonium and natural 
uranium (reactor according to para. b, Appendix) 


available burnup amounts to 1,32 n/kb and thus 1 ton of fuel will yield Q’ = 
= 6160 MWd heat energy. Excess heat energy — as compared to the preceding 
case — now amounts to AQ = 2810 MWd per ton. 

After the burnup level had been attained, the remaining Pu?®® quantity 
in the reactor started-up with plutonium fuel mixture will amount to 3,32 kg 
per ton of fuel. This is more, than the quantity of plutonium placed into the 
reactor for the start-up (2,2 kg Pu®*) and therefore, in the given reactor equi- 
librium fuel system appears to be realized, even when taking all possible losses 
into account. 

In a way similar to that given in the previous paragraph, use of relation 
(2) will define the value of plutonium, summarized by the following table : 


Preprocessing costs (pz) $/kg U3 25 SH 
Plutonium value (pp,) $/kg 11650 6660 1670 


Actual data are shown as a function of reprocessing costs by curve b 
of Fig. 6, proving that for the given fuel cycle and with the reactor as specified, 
a very considerable plutonium value is obtained, even in case of higher repro- 
cessing costs. This will indicate the optimistic attitude of interested experts 
towards the problem. (In this case the circumstance that a longer interval 


between plutonium production and utilization may unfavourably effect the 
valuation, has been ignored — see para. 5/d.) 
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5. Conclusions 


In order to compare the method of valuating plutonium, as proposed 
by the authors and based on energy production, against the method of 
valuating plutonium by its fixed price as known from literature, it seems 
to be advisable to introduce further refinements to the latter, rather arbi- 
trary procedure. 

As already mentioned in Section 1, plutonium as a concentrated fissile 
material produced in a reactor is sometimes evaluated in the literature as 
able to substitue the fissile material content of the U?®> — enriched uranium 
fuel (para. 1/a). For this evaluation, — specific data generally accepted by the 
literature will be assumed, according to which the cost of uranium enriched 
to 90 per cent U?> amounts to 15 000 or 30 000 $/kg according to enrichment 
costs [3]. However, for a more precise evaluation — if it is proposed to deter- 
mine the price of plutonium similarly with reference to the burnt-up material 
in the reactor — a few other factors must be taken into account, these latter 
being, as a rule, ignored by literature sources. These may be summarized in 
the following paragraphs : 

a) The thermal equivalents of U° und Pu®® or Pu*4! are by no means 
identical (see para. c, of the Appendix). 

b) The composition of the quantity (gp,) of plutonium produced, varies 
with irradiation or burnup level. This composition can be characterized by the 
isotope ratio (m) representing the ratio of fissile Pu-isotopes to the total quan- 
tity of Pu-isotopes. This is according to the analysis given in Section 3.1. a 
function of reactor type and irradiation level : 


ris iene) 7 TOO), (20) 
$(239+2407+ 241+242) 


c) Plutonium extraction costs (pp $/kg) must also be considered when 
valuating plutonium. 

Taking into account all the enumerated factors, the value of plutonium 
in the producing reactor fuel is given by 


tee pot =. Ske. (21) 
Spu 


In the above relation pp, denotes the price converted from the basic price 
(in $/kg) of enriched uranium in accordance with composition (m) and equi- 
valent heat of plutonium. Fig. 9 shows the plutonium values as computed 
from equation (21) as a function of irradiation. In this diagram plutonium 
values (pp,) are shown for the different U2%5 basic prices and different reprocess- 
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ing cost factors (pr). The diagram also contains a plot of the plutonium yield 


(gpy) and isotope ratio (m) against irradiation level. 
d) The value of plutonium is influenced by the time interval between 


production and utilization (discount value). The later plutonium will be utiliz- 
ed, the less its value is at the time of production. Assuming that plutonium 
extraction immediately precedes utilization, i.e. that there is no interest 
charged after pr, and if n denotes the number of years between plutonium 
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Fig. 9. Plutonium yield (gpu)s isotope ratio (m) and the value of plutonium calculated using 
fixed basic prices (Ppu) as a function of irradiation level 


production and plutonium utilization and p denotes the interest rate : then the 
value of plutonium in the burnt-up fuel of the producing reactor will be given by 


eg P u 
Pein oa Fh $/kg. (22) 


Fig. ° shows, for the sake of comparison, curves of plutonium valuation 
using the fixed basic price, but corrected according to the principles discussed 
(curve c, refers to a basic price of 15 000 $/kg, curve c, to 30 000 $/kg.) The 


fo ; ; 4 
ormer price, as seen from para. c, Section 1 is nearly consistent with the 
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price obtained on the fuel oil equivalence principle. [Burnup has been assumed 
identical with the cases given in Section 4 (Q = 3350 MWd/ton)]. As can be 
seen, the valuation curve according to the higher basic cannot be motivated 
by power engineering arguments, giving an unnecessarily high price to 
plutonium. Contrasts will become even sharper, if plutonium is produced by a 
different reactor, having e.g. a different initial conversion factor value or an 
initial excess reactivity, different from that assumed in our example. An ener- 
getically proper valuation will in this case show an even greater reduced pluto- 
nium value compared to the results obtained in our example. It is therefore, 
according to the opinion of the authors, a sound principle to avoid valuating 
plutonium — as a nuclear power cost factor — at fixed prices, taken arbitrarily 
or on the basis of any other considerations. Proper valuation must take into 
account design and operational characteristics of the plutonium-producing 
reactor, as well as the fuel cycle utilizing plutonium. It may be obvious that a 
similar energetic valuation possesses a lot of problems which should in the future 
be cleared up. Perhaps most important among these may be formulated in 
such a way that a solution most favourable for national economy should 
be sought for, by comparing different reactor types and fuel cycles. The proce- 
dure to be followed is similar to that, when during design of a conventional 
back-pressure type heat power plant the analysis is extended beyond the scopes 
of an industrial power plant and lifted to the level of national economy consid- 
erations [10]. Obviously, this will require more serious studies, but undoubt- 
edly the solution of these problems will advance the question of peaceful 
utilization of nuclear energy. 


Appendix 


a) Nuclear constants used for numerical solutions of the equations 
defining nuclear fuel composition variations (normalized data). 


SEE EE SEER 


ingh 4 Number of 

Thermal capture Thermal fission fission 

Isotope cross-section cross-section neutrons 
Gq barn oy barn 

1285 687 580 2,08 
U238 2,75 = is 
Pu239 1065 750 2,03 
Pu240 550 a ie 
Pu24t 1480 110 2,39 
Put? 30 — : 


ye 
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b) Characteristics of the hypothetical reactor used for the present analy- 
sis ; in hot start-up conditions : 


fuel natural uranium 


fast fission factor e = 1,030 
resonance escape probability p = 90,8828 
thermal utilization factor f = 0,865 
7 = 1,335 
Leakage factor P=(972 
initial conversion factor C, = 0.8 


initial excess reactivity 
for burnup and poisoning (with- 
out Xe and Sm) Oo; = 0,02 


c) The heat energy released during the fission of one gramme of fissile 
material : 


Isotope Energy MWd/ton 

(UE 0,94 

Pu? 0,878 

Put 0.79 
Summary 


The paper aims to consider one of the most discussed factors of the costs of electric 
energy produced in atomic power plants, i.e. to consider the value of plutonium produced. 
The basic consideration is that plutonium produced in an atomic power reactor is worth as 
much as the electric energy generated by it. This procedure is basically different from other 
means of evaluation — more or less arbitrary — discussed so far in the special literature, but 
it renders necessary to take into consideration when evaluating plutonium, the way of produc- 
tion in the plutonium producing reactor as well as the fuel cycle of plutonium utilizing. This 
procedure considers even such factors as arise from an eventually late utilization of the pro- 
duced plutonium. Thus it can be supplemented at any time — when more exact data are 
available — with expense factors arising from investment costs. To illustrate the course of 
calculations there is a more detailed examination of the burnup with prolonged fuel cycle, 
and further of the utilization of produced plutonium mixed with natural uranium in self- 
surtained cycle. Whereas the first procedure in view of the costs of plutonium extraction as 
known today, seems by no means to be economic, the economic conditions of the second 
procedure — under certain circumstances — are eyen nowadays available. The procedures 


published in the paper may serve as a basis for further discussions for a further solution of 
the problem. 
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EINIGE BEMERKUNGEN ZU EINEM KAPITEL 
DER THEORIE DER GITTERGESTEUERTEN GLEICHRICHTER 
Von 
O. BENEDIKT und L. ViTALyos 
Lehrstuhl fiir elektrische Spezialmaschinen an der Technischen 
Universitat, Budapest 


(Eingegangen am 15. August 1958) 


Ein nicht unwichtiges Kapitel der Theorie der gittergesteuerten Gleich- 
richter befaBt sich mit der Frage, wie sich der zeitliche Verlauf des gleich- 
gerichteten Stromes gestaltet, wenn der Gleichstromverbraucher nicht rein 
ohmisch ist, sondern auch Induktivitat und eine eigene elektromotorische 
Kraft enthalt. Kapazitive Lasten kommen im allgemeinen nicht vor, es darf 
also ein Verbraucher mit den geschilderten Eigenschaften als der allgemeinste 
Fall betrachtet werden. 

Der gleichgerichtete Strom wird durch die Differentialgleichung 


w—Uy= RTL +U, (1) 


bestimmt, wo wu die gleichgerichtete Spannung ist, d. h. jener Teil der Sekun- 
darspannung des den Gleichrichter speisenden Transformators, der durch den 
Gleichrichter an die Klemmen der Last geschaltet wird. R ist der Ohmsche 
Widerstand, L die Induktivitat, U; die innere Spannung der Last, wihrend U, 
die Bogenspannung bezeichnet, von der angenommen wird, daB sie von der 
GréBe des Stromes i unabhangig konstant bleibt. 

Im stationdren Betriebszustand ist die Spannung u periodisch, die Dauer 


2 7 
einer Periode entspricht dem elektrischen Winkel—— , wennm die Phasenzahl 
m 


des Gleichrichters bezeichnet. Bekanntlich kann der Gleichrichter in zweierlei 
Betriebszustanden arbeiten [1] : 


2% 
A) Wahrend der Gleichrichtungsperiode—— wird die Phasenwicklung 
m 


des Speisetransformators durch den Gleichrichter ununterbrochen an den 
Klemmen der Last gehalten. Da dies nur dann geschehen kann, wenn der 
Gleichrichter in dieser Zeitspanne auch Strom liefert, wird die geschilderte 
Betriebsart als der Zustand des kontinuierlichen (die aneinanderfolgenden 


2m ' 
—— Zeitspannen ununterbrochen ausfillenden) Leitens bezeichnet. 
m @ 
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B) Wird die Phasenwicklung des Speisetransformators durch den 
2% c 
Gleichrichter nur wahrend eines Teiles der Gleichrichtungsperiode an die 


Last angeschlossen, und erhalt der Verbraucher wahrend des iibrigen Teiles 
dieser Zeitspanne keine Speisung vom Gleichrichter, so ist die Leitung 
lickend. 

Ob der Gleichrichter in dem einen oder in dem anderen Betriebszustand 
arbeitet, bestimmen die Kenndaten der Last sowie der durch die Gittersteue- 
rung festgesetzte Ziindwinkel gemeinsam. 

In Bild 1 ist als Beispiel der charakteristische Verlauf der gleichgerich- 
teten Spannung u bei liickender und bei kontinuierlicher Leitung fiir den 


Bild 1. Verlauf der gleichgerichteten Spannung u bei aussetzender und bei nicht aussetzender 
Leitung 


Fall einer dreiphasigen Gleichrichtung (m = 3) aufgetragen. Die Kurven 
zeigen den allgemeinen Fall, in welchem die Ausgangsspannung des Gleich- 
richters im Verlauf einer Gleichrichtungsperiode (infolge der Induktivitat der 
Last) auch negative Werte annimmt. Im Bild wie auch in der folgenden Bespre- 
chung wurde die wegen der verzégerten Kommutierung entstehende Anderung 
der Spannungsform vernachlassigt, d.h. es wurde angenommen, daf die 
Kommutierung momentan erfolgt. 

Auch aus Bild 1 ist ersichtlich, da die Schwierigkeit der Lésung der 
Differentialgleichung (1), die auch in der Form 


Rit Lb = u—(U, +0) (2) 


geschrieben werden kann, darin besteht, daB die aus sinusférmigen Teilen 
zusammengesetzte Spannung u blo in je einem beschrankten Gebiet einer 


22 
Dauer von héchstens —— in geschlossener Form aufgeschrieben werden kann 
mw ; 


Dementsprechend hat auch der so erhaltene analytische Ausdruck des Stromes 
nur in diesen Gebieten Giltigkeit. 


— 
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Setzt man den Ursprung des Koordinatensystems in den Augenblick 
der Ziindung der einen Anode (oder der einen Réhre), d. h. in den Anfangspunkt 
einer Gleichrichtungsperiode, und rechnet man den Ziindwinkel bei beliebiger 
Phasenzahl vom Beginn der Speisespannungs-Halbperiode, laBt sich Gleichung 
(2) in folgender Form schreiben : 


Ri+Lo = U,, sin (wt + a,) — (U; + Us). (3) 


Nimmt man an, daB die Induktivitat L vom Strom unabhangig ist, 
wird die Differentialgleichung eine lineare mit konstanten Koeffizienten. Es 
soll ferner die Annahme gelten, da die innere Spannung des Verbrauchers 
U, konstant ist. 

Wie aus der Theorie der linearen Differentialgleichungen mit konstanten 
Koeffizienten bekannt, besteht die Lésung aus zwei Teilen, aus dem stationa- 
ren Strom 
U,+U, 


= (4) 


i, = —B sin (ot + a, — 9) te 


und aus dem transienten Strom 
(5). 


Mit den neu eingefiihrten Bezeichnungen 
P= ier und Z= R? + @L) 
schreibt sich also die vollstandige Lésung in der Form 


R 
U; + Uo Fis (6). 


i= i,t i, =H sin (01+ a, — 9) — = + ke 


Die Gleichung wird folgendermaBen umgestaltet : 


1 
statt — setzen wir— COS, 
Z R 


Un 
fa ausgeklammert, 
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U,+U as 
die pociiene = a eingefuhrt, 
Us 
R ; at 

und statt des Exponenten — L t des Exponentialgliedes setzt man ee ’ 
worauf man 

mens ; =e ; 

edge CLA Aa NCA Be ] (7) 


erhalt. 


Die fiir die liickende und fiir die nichtliickende Leitung giiltigen 
Lésungen der Gleichung (7) unterscheiden sich in den Anfangsbedingungen. 
Bei liickender Leitung ist der Stromkreis im Moment der Ziindung stromlos, 
es ist somit i = 0, wenn wt = 0. Wird dieser Wert in die Gleichung (7) ein- 
gesetzt, erhalt man 


0 = cos@ sin (a, — y) —a+ K, (8) 


woraus 


K = a — cos 7 sin (a, — @). (9) 


Setzt man den so erhaltenen Wert von K wieder in die Gleichung (7) ein, dann 
wird 


wt 


en U te 
te i {cos @ sin (wt + a, — y) —a-+ [a — cos @ sin (a, — ¢)]e tor) . (10) 


Die Gleichung beschreibt den Stromverlauf richtig, solange der daraus berech- 
nete Strom t > 0. Da der Strom nur in einer Richtung flieBen kann, wird 
er, nachdem er bei dem Phasenwinkel wt = a, verschwunden ist, bis zur Ziin- 
dung der nachsten Anode (Réhre) den Wert Null beibehalten. 

Bei nichtliickender Leitung mu8 die Lésung der Gleichung (7) die 


Bedingung befriedigen, da®B der Strom im Zeit nine nach der Ziéindung 
mw 


denselben Wert annihme, wie im Zeitpunkt der Zindung, da in stationiarem 
Zustand und bei symmetrischer Schaltung die in der Leitung an die Reihe 
kommende Anode (Réhre) von der vorangehenden Anode den gleichen Strom 
ubernimmt, den sie der nachfolgenden ubergibt. Bezeichnet man diesen Strom 


mit I,, so nimmt diese Bedingung unter Anwendung der Gleichung (7) fol- 
gende Form an 
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Ure ; 
Tee Bicosy ale apenas cos @ sin eas — 
R R m 


20 


ase ke: (11) 


Nach Ordnen wird hieraus 


sin 


22 } 
“= + a9) — sin (a. —#) 


m 
splat id aa ’ (12) 


was mit Hilfe der aus der Trigonometrie bekannten Beziehungen auch in 
folgender einfacherer Form geschrieben werden kann: — 


2 cos p 


2a 
eee oe 


K = 


sin —- cos [ —gpt = ; (13) 
m m 


Setzt man den so erhaltenen Wert von K in die Gleichung (7) ein, bekommt 
man folgenden Ausdruck fiir den Strom : 


c= [erinert ans + 


wt 


as ecard aie a 
m 


Fee CH 
sin cos 
2% m 


pat eee 


der fir den Fall der nichtliickenden Leitung giiltig ist. 

Es muB nun noch jene Grenze festgestellt werden, bis zu der die Glei- 
chung (10), und ab welcher die Gleichung (14) giiltig ist, bei der also nach dem 
Betriebszustand der liickenden Leitung der Betriebszustand der nicht- 
liickenden Leitung beginnt. Diese Grenze liegt dort, wo bei der Ziindung 
der in der Leitung an die Reihe kommenden Anode (Réhre) die vorangehende 
Anode (Réhre) eben nicht mehr stromlos ist, oder wo die Brenndauer der 


2% 
Anode (Réhre) genau—— bhetragt, d.h. eine volle Gleichrichtungsperiode 
mw 


umfaBt. 

Bei liickender Leitung ist die Dauer der Stromfihrung einer Anode 
(Réhre) in elektrischem Winkel gemessen a;—a,, und a, kann aus der Glei- 
chung (10) auf Grund der Bedingung berechnet werden, da® der Strom ver- 
schwindet : 
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aj[—az 


pene cosy sin (a, — a, + a, — p) —a + [a — cos psin (a, — p)je ?}. 
R 7 


(15) 
Nach Ordnen 1aB8t sich die Bedingung (15) auf folgende Form bringen: 


a; a, 


[cos @ sin (a, — y) — a] e'é? — [cos ysin (a, — y) —a] eee, (16) 


Die Lésungen a, # a; der transzendenten Gleichung geben den Léschwinkel 
a, als Funktion des Ziindwinkels a, und der Parameter gy und a an. 
Ein Nomogramm, das die nach a und cos y geordneten Lésungen der 
Gleichung (16) enthalt, wurde 1945 von K. P. Pucutowsxt veréffentlicht [2]- 
Der Betriebszustand der liickenden Leitung besteht, solange der aus 
der Gleichung (10) berechnete oder der aus dem erwahnten Nomogramm 


2% 3 : 
abgelesene Wert a, kleiner ist als 4, + . Erhalt man fir a, einen gréBeren 
m 


Wert, besteht bereits der Betriebszustand der nichtliickenden Leitung, 
und dementsprechend wird a, nicht mehr dem aus Gleichung (10) berechneten, 


27 
sondern dem Wert @ + — gleich sein. 
m 


Die soeben geschilderte, im Schrifttum allgemein tbliche Methode, die 
Gleichungen (10) und (14) abzuleiten, die den zeitlichen Verlauf des durch den 
gittergesteuerten Gleichrichter gelieferten Stromes in den beiden Betriebs- 
zustanden beschreiben, ist in gewissem Sinne unvollkommen, trotzdem sie 
sowohl vom physikalischen als auch vom mathematischen Gesichtspunkt aus 
als korrekt bezeichnet werden mu. Dieser gewisse Mangel besteht darin, 
dafs zwischen den beiden der Form nach voneinander stark abweichenden 
Gleichungen im Laufe der Ableitung keine Korrelation zustande kommt, 
die klar zeigen wiirde, wie die beiden Gleichungen an der Grenze der nicht- 
liickenden Leitung ineinander itibergehen. 

Von einem der Verfasser wurde eine Methode ausgearbeitet [3], die den 
erwahnten Mangel der beschriebenen Ableitung teilweise aufhebt, diese 
Methode wurde aber in erster Linie zur Berechnung der infolge Anderungen 
des Ziindwinkels auftretenden transienten Stréme ausgearbeitet ; sie bedient 
sich dementsprechend mathematischer Methoden, die bei der Behandlung 
transienter Erscheinungen zweckmaBig angewendet werden [4]. 

Der Gedankengang der Methode ist folgender: Die gleichgerichtete 
Spannung wu kann in beiden Betriebszustanden in Spannungsimpulse zerlegt 


Ut eae 2m 
werden, die sich in Zeitspannen T= -—— wiederholen. Jeder Impuls 1aBt 
mw 


sich in Koordinatensystemen, deren Nullpunkte der Reihe nach in den Zeit- 
punkt der Ziindung verlegt werden, der Gleichung u = U,, sin (wt + a,) 
gem4$ auftragen, u. zw. bei liickender Leitung im Bereich 0 < ot < [— Az» 
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bei nichtliickender Leitung hingegen im Bereich 9< wie . Hierauf 
hért der Spannungsimpuls auf und es folgt ein neuer. - 

Ist die Stromleitung liickend, entstehen die einzelnen Stromimpulse 
durch die Einwirkung der einzelnen Spannungsimpulse unabhangig von den 
vorangegangenen. Jene Spannungsimpulse, die sich an den Klemmen des 
Stromkreises schon abgespielt haben, wirken in der Erzeugung des Stromes 
nicht mehr mit, da der gespeiste Stromkreis infolge des Absinkens des Stromes 
auf Null nach jedem Impuls unterbrochen wird. 

Ist die Stromleitung nichtliickend, nimmt der Strom der einzelnen 


_ Anoden (Réhren) nicht den Wert Null an, bevor noch die in der Reihenfolge 


nichste Réhre ziindet. Man kann demzufolge die einzelnen Stromwellen nicht 
als in sich abgeschlossene transiente Vorgange berechnen, in denen die voran- 
gegangenen Spannungsimpulse auf die Form der Stromwelle keine Wirkung 
ausiben. 

Wahrend nimlich der Stromkreis, wie schon oben erwahnt, bei liicken- 
der Leitung infolge des Verléschens der Anode (Réhre) nach jeder Stromwelle 
unterbrochen wird, bleibt der Gleichstromkreis bei nichtliickender Leitung 
standig geschlossen. 

Fs bezeichne i* den transienten Strom, der in einem Stromkreis mit dem 
Ohmischen Widerstand R, der Induktivitat L und mit einer Stromquelle der 
inneren Spannung U; + U, entsteht, sofern man an den stromlosen Kreis die 
Wechselspannung u = U,sin (t+ a,) bei einem Phasenwinkel wt = 0 
anlegt. Diesen Strom beschreibt Gleichung (10), es ist somit 


Ot 


tee 


r= 2 cos sin(wt + a,— y) —a+ [a — cos psin (a, — ¢)]e : ( (17) 


Ist die Stromleitung liickend und entsteht jede Stromwelle unab- 
hangig von den vorangehenden, ist der Strom gleich i* im Bereich 0< ot < a;—a@,, 
wo a, aus der Gleichung i* (a, — az) = 0 ermittelt werden kann. a, ist selbst- 
verstandlich die Lésung der Gleichung (16). 

Ist die Stromleitung nichtliickend, zerlegt der Rechnungsvorgang 


2 7 ; : 
den analytischen Ausdruck einer—— langen Periode des gleichgerichteten 
m @ 


Stromes in zwei Teile : 

1. in den durch den soeben wirkenden Spannungsimpuls erzeugten 
Strom i; und 

2. in die Superposition der durch samtliche vorangegangene (in dem 
untersuchten stationaren Zustand unendlich viele) Spannungsimpulse erzeug- 
ten Stréme, die mit i;; bezeichnet werden sollen. 

Der Strom i; ist gleich i*, i;; 1aBt sich dagegen durch folgende Uber- 
legung ermitteln. 
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Die Stromkomponente, die der dem soeben sich abspielenden Spannungs- 
impuls vorangehende n + 1-te Impuls erzeugt, ist : 


yes (t-+nt) 
i Rye (18) 
woraus 
; co — a (t+n7) 
HESS NaC ery: (19) 


== 


Hebt man die von n unabhangigen Faktoren aus (19) heraus, erhalt man 


= t co nt 2 
Hier Sp (20) 
Da ries lee 1, 148t sich schreiben: 
CerE: 1 See 
yy = 1 (t) 3 ts (21) 
1 — oe? 


Fir den vollen Strom bei nichtliickender Leitung gilt somit 


: 4 3 : k 1 sei? 
t= i +i = + (7) eect Aes (22) 


== 
Theres 


Mit Verwendung des fiir i* erhaltenen Ausdruckes (17) kann man sich 
davon tiberzeugen, daB Gleichung (22) mit Gleichung (14) identisch ist. 

Wie man sieht, ist es dank der oben kurz geschilderten Methode gelungen, 
zwischen den Ausdriicken der liickenden und der nichtlickenden Leitung 
einen engen Zusammenhang herzustellen. (22) zeigt, daB bei nichtliickender 
Leitung der Ausdruck des durch den eben wirkenden Spannungsimpuls ent- 
standenen Stromes derselbe ist wie bei liickender Leitung, daB® ferner eine 
durch samtliche friihere Spannungsimpulse hervorgerufene, exponentiell 
abklingende Komponente diesen Strom derart modifiziert, daB die Symme- 
triebedingung i(0) = i(t) erfillt wird. Aus der Gleichung (22) folgt 

(0) =i@) =F) -—_, (23) 
Leereutet : 
da 


i* (0) =0. 


7 
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An der Grenze der liickenden und der nichtliickenden Leitung ist 
i*(t) = 0 und daher i =i’. 

Das hier beschriebene Berechnungsverfahren schafft eine viel iiber- 
 sichtlichere mathematische Korrelation zwischen den beiden Betriebszustanden 
als die klassische Methode der Lésung der Differentialgleichung. Der Gedan- 
 kengang der Ableitung leidet jedoch an dem Mangel, eine ziemlich kompli- 

zierte, auf dem Gebiete der Berechnung transienter Stréme itibliche Zer- 
legung zur Berechnung eines im wesentlichen stationaéren Vorganges zu ver- 
wenden. Der Vorzug der angefihrten Methode zeigt sich nur dann, wenn man 
nicht den stationaren Strom zu berechnen wiinscht, sondern voraussetzt, 
_daB die Stromleitung um endliche Gleichrichtungsperioden friiher begonnen, 
‘oder der Ziindwinkel sich um ebenfalls endliche Gleichrichtungsperioden 
friiher geandert hat. 

Um diese Mangel der bisher bekanntgewordenen Methoden zu beseiti- 
gen, wird im folgenden ein neues Verfahren vorgeschlagen und erlautert. 
Gema®R der Meinung der Verfasser schafft es einen einfacheren und physika- 
lisch anschaulicheren Zusammenhang zwischen den beiden Relationen (10) 
und (14). Diese Ableitungsmethode stimmt bis zur Ermittlung des Zusammen- 
hanges (10) mit der Klassischen Methode der Lésung der Differentialglei- 
chung iiberein, um von da ab abzuweichen und sich nicht mehr auf die un- 
mittelbare Bestimmung des Wertes K aus den Grenzbedingungen wt = 0, 
Cel, aOtS = ,i = I, zu stiitzen. Die Ableitung verlauft wie folgt : 

Es sei angenommen, daf die fir die nichtliickende Leitung giltige 
Stromform durch den vorlaufig noch unbekannten Zusammenhang t = 1(at) 
beschrieben wird. Diese Gleichung mu8 — wie bekannt — die Bedingung 
erfillen, da®8 der Strom am Anfang und am Ende der Leitungsperiode einer 
Anode (Réhre) den gleichen Wert J, annehmen soll. 

Die von den Verfassern vorgeschlagene Ableitung geht davon aus, daB 
der Zustand der liickenden Leitung aus dem Zustand der nichtliickenden 
Leitung nicht nur in der bisherigen Weise erreicht werden kann, da8 bei 
konstanter innerer Spannung U; der Ziindwinkel a, vergréBert wird, sondern 
auch dadurch, daB bei konstantem Ziandwinkel a, die innere Spannung U; des 
gespeisten Stromkreises zunimmt. 

VergréBert man den Wert U; bei konstantem Ziindwinkel a,, wird der 
durchflieBende Strom standig kleiner, da die Richtung von U; dem Stromflu8. 
entgegenwirkt. 

Der entstehende Strom kann auf Grund des Prinzipes der Superposition 
auch als die Differenz des durch die Spannung u(t) verursachten pulsierenden 
positiven Stromes und des durch U; + U, verursachten negativen kon- 
stanten Gleichstromes aufgefaBt werden. VergréBert man U;, so nimmt die 
negative Gleichstromkomponente zu. Die Wechselkomponente des Stromes. 
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wird durch die Wechselkomponente der an den Stromkreis angeschlossenen 
Spannung u(t) bestimmt. Da die Spannung u(wt) trotz der Anderung von 
U; unverandert bleibt, solange die Stromleitung tatsachlich nichtliickend 
ist, bleibt auch die Wechselkomponente des Stromes unverandert. Erreicht 
U; einen gewissen Wert U;, ist der Absolutwert der negativen Stromkompo- 
nente gerade um I, gréBer geworden, der Strom am Beginn und am Ende der 
Leitungsperiode einer Anode (Réhre) ist somit gleich Null. Dies ist offenbar 
die Grenze der nichtliickenden Leitung bei dem gegebenen Ziindwinkel 
a, und dem gegebenen Phasenverschiebungswinkel gy. Fir diesen Zustand 
erhalt man statt des friher definierten allgemeinen Wertes a einen der GréBe 
U; entsprechenden Wert 


jae seo (24) 


Ist der Strom hier an der Grenze der nichtliickenden Leitung durch 
die Gleichung i’ = i’(wt) gekennzeichnet, ist der der urspriinglichen inneren 
Spannung U; zugeordnete »nichtliickende« Strom 


i= +]y, (25) 


U 


Ue U,) ab hn 
eae <9 Ee Be) = (a — 0) =", (26) 


R 


da der Gleichstromwiderstand des Stromkreises R ist. Der Strom i’ kann, da 


es sich hier um einen Grenzwert handelt, noch aus dem Zusammenhang (10) 
berechnet werden, wenn man a’ statt a setzt. 


Somit ist zwischen der fiir die liickende Leitung giiltigen Gleichung (10) 
und der den nichtliickenden Strom beschreibenden Gleichung (18) ein enger 
und anschaulicher Zusammenhang geschaffen, ist es doch ersichtlich, daB der 
Grenzfall i’ des durch die Gleichung (10) beschriebenen Stromes einen Teil des 
durch Gleichung (14) angegebenen Stromes darstellt. 


Z, oo . 2 of oo, ? 
sass ur Vatsachienen Berechnung vonz und I, benétigt man den Wert von 
a dessen Ermittlung Gleichung (16) erméglicht, wenn man bericksichtigt, 


daB an de r ul i 28 
r Grenze der nichtliickenden Leitung @ = 4,+ —, undsomit 


20 
a 
2 m az 
ter 


4 22 Ls 
Joos psin «+= — 9] ~a'le pe = [cos p sin (a, — p) —a'Je : oo, (at) 


m 


é ne mili 


se 
7 
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oP 


Dividiert man beide Seiten durch ets? und lést man die Gleichung auf a’, 
erhalt man: — 


sin (a, — y) —sin}a, + cable j ene 
’ m 
a => cos P a (28) 
al pee gue 


Nunmehr 1a8t sich auf Grund der Zusammenhange (10), (25), (26) und (28) 
der Stromverlauf fiir jeden beliebigen konkreten Fall aufschreiben. Dazu ist 
nur folgendes nétig : man berechnet a’ aus (28), und daraus gem48 (26) den 


Bild 2. Verlauf des Stromes, der in einem Stromkreis mit einer Zeitkonstante L/R = 0,01 sec, 
bei m = 3, ay = 70° und bei verschiedenen Werten von a flieBt 


Wert I, ; auBerdem setzt man in Gleichung (10) statt a den soeben erhaltenen 
Wert a’ ein und bekommt derart den Ausdruck des Stromes i’. Addiert man 
gemaB Gleichung (25) i’ und I,, erhalt man den vollstaéndigen Ausdruck fur 


den nichtliickenden Strom. 

In Bild 2 ist beispielsweise der zeitliche Verlauf des Stromes dargestellt, 
der in einem Stromkreis mit der Zeitkonstante L/R = 0,01 sec, bei m = 3, 
a, = 70°, und bei verschiedenen Werten von a flieBt. Fir a’ ergibt sich aus 


Gleichung (28) fiir diesen Fall ein Wert von 0,58. 


Anhang 


Setzt man die Ausdriicke fiir i’ und I, aus den Gleichungen (10) und (26) 
in Gleichung (25) ein, und verwendet man auch den Ausdruck (28) fiir a’, 
mu man natirlich Gleichung (14) erhalten. 


3 Periodica Polytechnica El II/4. 
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Zur Kontrolle soll diese Berechnung durchgefiihrt werden : 


j = {eos psin (ot +a, —9)— a + [a’ — cos psin (a, — y)Je*”} + 


+52 (a —a) = 


ot 


= A feos y sin(wt+a,—) —a-+ [a’ — cos@ sin (a, — ¢)] e. ee - (29) 


Mit dem Ausdruck (28) fiir a’ in Gleichung (29) wird der Koeffizient 
des exponentiellen Gliedes : 


mtg @? 


: 2 
sin (a, — ) — sin (a, += — rle 
m 


—sin(a,—g)}. (30) 


2% 


leu 


Wenn man auf gemeinsamen Nenner bringt und zusammenzieht, liest man 


sin (a, — y) —sin|a, MS o| lems" 
A = R cos Pp aes (31) 
1—e"= 


22 
Dividiert man Zahler und Nenner durch —e™‘*? und verwendet man den 
trigonometrischen Satz von der Umformung der Differenz der Sinusse zweier 


Winkel in ein Produkt, so ergibt sich folgender Ausdruck : 


oe 2 cos p 
R 2 


‘ee ay mtgp 


eee 
sin cos 
m 


e+). (32) 
m 


Nach Einsetzen dieses Ausdrucks (32) in (29) erhalt man tatsichlich den Zusam- 
menhang (14): 


i Ops : 
i= P| cospsin(we+ a, — 9) — a+ 
2 cos p . eee 
+ BSE an Zcole, 4 2] i 
l—» mee a 
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Zusammenfassung 


Es werden zunichst diejenigen Methoden zusammengefafit, die man im allgemeinen 
zur Ableitung der mathematischen Ausdriicke fiir den von einem gittergesteuerten Gleichrich- 
ter gelieferten Strom verwendet, wenn der Verbraucher nicht rein Ohmisch ist, sondern auch 
Induktivitat und eine eigene elektromotorische Kraft enthalt. Est wird der Mangel dieser 
Berechnungsmethoden betont, die Betriebszustande der lickenden und der nichtliik- 
kenden Leitung voneinander getrennt zu behandeln und zwischen den Gleichungen, die 
diese beiden ineinander iibergehenden Betriebszustinde kennzeichnen, keinen anschaulichen 
Zusammenhang zu schaffen, oder diesen nur durch ein ziemlich kompliziertes Verfahren mit 
betrichtlichem mathematischem Aufwand zustande zu bringen. Die Verfasser beschreiben 
ein einfaches Verfahren zur Herstellung einer engen und physikalisch sehr anschaulichen 
Korrelation zwischen den die beiden Betriebszustande kennzeichnenden mathematischen 
Ausdriicken. 
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I. Introduction 


The present stage of development of automatic telephone equipment is 
marked by the headway the use of crossbar switches is making. Unlike friction- 
brush type switching equipment, with brush assemblies sliding along a straight- 
line or curved tracks, modern crossbar switches operate only with relay-like 
motions. Consequently, a substantially lower amount of power is required for 
the operation of this new-type switch than is needed for earlier types of equip- 
ment. Obviously, the question may be raised as to the optimum value of the 
operating voltage of telephone exchanges built up on crossbar switches. 

This paper examines the reasons why the voltage of batteries used in 
manually operated telephone exchanges had to be increased from the original 
24 volts to 48, or even 60 volts, on the introduction of automatic service. 
A survey will also be given of major considerations determining the rating 
of voltage. 

As is known magnets and relays account for about sixty per cent of 
the overall costs of a telephone exchange. It is essential, therefore, that by 
considered assessment of the different factors a voltage value should be deter- 
mined for the operation of a particular exchange, which in turn might be apt 
in resulting the possibly highest reduction of costs. 


Il. Calculation of windings 


For reasons of economy relays and functional magnet assemblies should, 
in general, be of a size that a minimum quantity of copper should suffice for the 
generation of the required pulling force. Exception may, however, be made 
for cases when conditions bearing on the duration of operation, or such as are 
relate to a pre-determined resistance value, have been specified for the partic- 
ular magnet assembly. 

The operating force may be written as 


Poe 10 bt 
2uF 
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where F denotes the cross-section of iron in sq. centimetres, 
the magnetic permeability of iron, 
® the magnetic flux obeying the equation 


Db — BEIOEE 10-8 volt sec 
ae 0 
LE, Mok 


where ni is the number of field ampere turns, 
1, the length of the iron portion of the magnetic circle in centimetres, 
|, the length of the spark gap in centimetres, 
F,, the cross-section of the air gap in sq. centimetres, a value as a rule 
larger than that of F;. 
For a given type of magnet, i.e. for specified geometrical dimensions, 
@ is the function of ni, while the value of P depends of that of (ni)?. Conse- 
quently, when the winding space has been defined, a maximum number of 
ampere turns may be obtained by the least amount of copper, provided that 
the copper had been loaded with a maximum current density. 
The current intensity obtainable in the winding may be written as 
follows : 


eOlTT 
Jas 


10? amps 


where 1, denotes the maximum permissible current intensity in the winding 
expressed in terms of amps/sq. millimetre, for windings generally 
employed in telephone technique, and at a continuous load of the 
order of 4 to 5 amps/sq. millimetre, 
6 the diameter of the wire in centimetres. 
Accordingly, the number of turns on a fully wound core, may be expressed as 


4Q 


62 2 


n=k 


where Q is the cross-section of the winding space in sq. centimetres, 
k a space filling factor, being the function of 6; for wire diameters 
current in telephone technique, i.e. diameters ranging from appr. 
0,08 to 0,4 millimetre, this factor, with increases in the value of é, 
may become 0,45 to 0,6. 


The number of ampere turns obtainable may be expressed by the follow- 
ing equation : 


: 4Q . &x ’ 
nt oe Re aD ees 
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Since both @ and i, are constants, the equation is reduced to 
ni = const k 


i. e. for a definite winding space the maximum number of obtainable ampere 
turns is independent of variations of the two components of ni, i. e. the value of 
ni only increases when the value of k is sufficiently large, or, with other words, 


mA Vin 
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Fig. 1. The parameters of relays as functions of the diameter of coil wires 


a larger number of ampere turns can be obtained for windings made of sufficiently 


thick wire, i.e. windings having a lower resistance value. 
The value of the voltage required for the specified current intensity 


is expressed by the following equation : 


E=i.-R, volt 


where R, is the resistance of the winding expressed in ohms, t. e. 
Bee, 
4 8x dx 


R,=1, 
Im being the average length of a single turn on a fully wound core 


E = const: i volt 
82 


I. e. while the increase of the current intensity can be expressed as the 
meter, the voltage drops by the inverse square 


square ratio of the wire dia 
This means that the power consumption increases 


ratio of the wire diameter. 
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with the wire diameter only insofar as the value of k increases. Consequently, 
the curves of both ni and the power consumption are essentially identical with 
that plotted for k as a function of 6, as shown in Fig. 1. 

For purposes of practical calculations certain limitations may enter the 
picture. Thus the current intensity possibly should not exceed a value of 0,3 
to 0,5 amp. lest the point-like contacts of precious metal should rapidly deterio- 
rate under the impact of the contact pressure of 20 to 30 grammes in general 
characteristic of relays, particularly at breaking inductively loaded circuits. 
At a continuous load the watt consumption of relays should not rise above 
3.5 to 4 watts, lest the accumulated heat should have a damaging effect on 
the insulation of the enamelled wires. 

There is a minimum size for the construction of relays, determined by 
the length of springs employed. Furthermore, the volume of the relay will 
depend on the maximum number of contacts required, and on the number of 
ampere turns specified for their operation. The maximum number of contacts 
are the function of the circuit requirements. The so-called flat type R relays in 
current use, in this country, may have as many as fourteen contact springs 
(e. g. two change-overs and four makes), the cross-section of the winding space 
being approximately three sq. centimetres, the volume of the winding about 
12 cubic centimetres. With a relay winding of this size, and a wire diameter 
of 0,1 millimetre 


ni = k-Q-i, = 0,45 - 3,500 = 675 


ampere turns can be obtained, while with a thicker wire of, let us say, 0,4 milli- 
metre in diameter, the number of ampere turns can be increased to about 900. 
The number of ampere turns obtainable for each cubic centimetre of 
winding is of the order of 55 to 75. 
Calculations show that in the previous case the resistance of the relay 
is 1650 ohms, the number of turns being 18 000. The voltage required for oper- 
ation may be written as 


B= Re 65h 
Tt 


The resistance of a winding made of a 0,4 millimetre wire is 8 ohms, the 
number of turns being 1400, so that for 900 ampere turns a voltage of 


900 
E> = 8 = 5.2 volts 
1400 
will have to be specified. 


In this latter case current intensity takes on a too large value. If the 
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current intensity is limited to e. g. 0,1 amp, then in order to ensure a current 
density of 5 amps per sq. millimetre the wire diameter will have to be 


6 = 0,16 millimetre 


With a wire of this diameter the resistance of the realizable winding will be 
285 ohms, and the number of turns 7800. The maximum number of ampere 
turns obtainable reads 

ni = 780 
and the voltage required 


1 180 985 = 28 volts 
7800 


(The consumption of the relay is 2,8 watts.) 

For all practical purposes, with now current designs 200 to 250 ampere 
turns generally suffice for the operation of the fourteen contact springs, men- 
tioned earlier in this section. When no operating time has been specified, with 
a reasonable margin of safety 300 ampere turns have to be provided for the 
operation of the relay, so that a voltage of less than 28 volts will suffice. 

As a matter of course relays used in more complex circuits may operate 
in a number of combinations, and therefore they may have to be provided with 
two, or even more, windings. On the other hand, in the majority of cases the 
number of contact springs to be provided for relays is substantially below 
fourteen. However, in equipment using crossbar switches at several places, 
e. g. for the connection of the common marker circuits, relays with a large num- 
ber of make contacts may be required. The design of these types of relays 
can by no means be identical with that of the active relays of the circuits, 
although the magnetic pull force in this type of relays may on hand of the 
formula derived earlier in this action be even 1000 grammes, provided that 
soft-iron of satisfactory quality has been used. This force is appropriate for 
the operation of twenty to twenty five make contacts. 

The battery of the exchange will have to be designed for the supply of 
a voltage at which the relays may in a satisfactory manner assist in the solution 
of the overwhelming majority of circuit problems. Obviously, in particular cases 
two relays may have to be employed for the solution of certain specific prob- 
lems. For reasons of economy the design of the relay will have to be adapted 
to conditions most frequently encountered. 

The preceding calculations tend to confirm that even with current types 
of relays the number of ampere turns most frequently satisfying recurring 
requirements may be accomplished with voltages below 48 volts. It was not due 
to relays that in automatic exchanges the voltage had to be raised to values 
above 24 volts, the usual value for manual equipment. 
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In practice for a large percentage of relays the number of contacts 
may be less than that of the maximum carrying capacity specified for the 
particular relay. Consequently, a lower number of ampere turns will suffice 
and with the customary thin wires the resistance needed will not be ensured 
at a given voltage. In this case part of the winding may have to be replaced 
by resistance wire. The larger the battery voltage than the voltage value at which 
the majority of relays is realizable without the substitution of resistance wire, 
the larger the number of relays for which at an increasing rate resistance wire 
will have to be used. Obviously, this condition will add to the costs of relays. 

Stepping or switching magnets may occasionally have to exercise a 
substantially larger pull force than relays. In this case a larger winding volume 
is required to ensure an appropriate number of ampere turns. In the Strowger 
system, and other similar systems operating on the step-by-step principle, 
the size of the windings is determined by requirements set by the operating 
time rather than by other conditions. In systems operating with crossbar 
switches the magnets will have to overcome the resistance of the springs, and 
the pull force to accomplish this does not in general exceed the values obtain- 
able with windings of a size of those of the relay type. However, here the time 
factor may enter into the picture, and substantiate the use of larger voltage. 


III. Operation time 


The electrical portion of the actuating time of a relay, i. e. the time laps- 
ing before the armature starts, may be written as 


t= —' In sec 


where L,/R, is the time constant of the winding, 


I, the setting current intensity at which the armature starts, 
I, the definitive current intensity resulting from the ohmic law. 
The value of the time constant may be expressed as 


0,4 7 n2 


f= 


-10-8 henry 
h 0 


“,F, My Fy 


and for a given design as 


L, = const -n? henry 


4-k-Q\2 (6% x)2 k-Q 
- = const 
622 


16-k-Q-o-l,, 1,0 


L,/R, = const 


= const k 
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This means that time constant L,/R, is the function only of the space 
filling factor k. Essentially the actuating time depends on the second part of 
the equation, viz. 


The value of this formula is a function of the relation of I,,/I,. If 


I,/Ifg ——> 0 then t=0 
and 


T,/I[g — 1 then t=oco 


In the formula of the actuating time the expression In is essentially 
a multiplier. In practice for reasons of safe functioning 


ds maxe 058.45 


On the other hand the minimum is determined by a ratio of the minimum 
number of ampere turns required for the operation of a single make-contact. 


05 1 Im/lo 
Fig. 2. Operating time increase of relays as a function of their adjustment current intensity 


and the maximum possible number of ampere turns. The 70 to 80 ampere 
turns required for a make-contact represent approximately ten per cent of 
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the maximum number of turns obtainable. In practice, therefore, the value 
of the multiplier may vary between 
1 
. gh Reawres 
and 1— 0.8 
1 
In Se eee 
1—0.1 


1,6 
0,1 


This variation expressed as a function of I,,/I, is shown in Fig. 2. The value of 
the multiplier will become the unity, 1. e. the time constant will remain the 
sole factor influencing the actuating time, when I,, becomes about 63 per cent 
of I). 

Inasmuch as for a given winding space time constant L,/R, essentially 
varies with the space filling factor only, and is lower for thin wires, the only 
feasible scheme for shortening the actuating time is to increase the value of I,, 
I,, being of a definite value owing to the function to be performed by the mag- 
net. This is the reason that entailed the necessity of increasing the voltage. It was 
on this principle that the windings of the Strowger system of 60 volts 60 ohms, 
causing the brush assemblies of the selectors to step with a power consump- 
tion of 60 watts, had been developed. Owing to the known properties of the 
Strowger system it was essential to ensure high-speed operation of the equip- 
ment. (Direct actuation in response to the dial pulses, free selection of limited 
duration between two digits, etc.). Magnets of this type are encumbered with 
a number of disadvantages, e. g. fire risk, interruptions of the current of one 
ampere equalling to heavy strains on the contacts, strong shocks causing noises 
towards the battery, complications involved in spark quenching, etc. 

Owing to the use of registers, conditions in the Rotary system are by far 
more favourable. When a 48 volt battery is used, 10 watt windings amply 
suffice for the operation of switching magnets. Improvements in the operating 
conditions of exchange equipment using crossbar switches are still more 
remarkable. Here time-conditioned, high-speed magnets can be dispensed 
with. The number of ampere turns required for operation may be accomplished 
with voltages below 48 volts. 

Since the size of relays, and in case of crossbar exchanges, of the magnets 
actuating the bars, does not even require 48 volts, less the 60 volts, certain 


other considerations will have to be examined which may have a bearing on 
the battery voltage. 


IV. Other considerations 


1. Problems of the feeding bridge 


In order to ensure an adequate microphone current, even for long lines, 
the ohmic resistance of the feeding circuit has to be kept as 


low as possible, 
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or in other words, the resistance of the feeding bridge determines the maximum 
permissible resistance of a subscriber’s loop. The impedance of the feeding 
bridge has to be large enough to check the speech currents. Furthermore the 
feeding bridge has to be free of fire risk, even under the most rigorous condi- 
tions. 

Fig. 3 shows the standard feeding bridge of the 7A2 Rotary equipment. 
The high impedance is ensured by a separate choke coil, the signalling relay 
being of the usual flat type. The figure shows the values of the ohmic resist- 
ance for a 48-volt battery. Calculations indicate that pure earth leakage on 
wire “b’”’ will cause a consumption of 5,5 watts only in the relay. 

For microphones now in use a minimum current intensity of 0,030 amp. 
is required. The upper limit is of no particular interest, partly because even 


Fig. 3. 7 A2 type feeding bridge 


for a 0-ohm loop the intensity will fail to reach the dangerous value of 0,1 amp., 
partly because modern telephone apparatus are provided with current limit- 
ing features. 

The microphone current of 0,030 amp. permits the following loop resist- 
ance maxima, provided that the resistance of the supply relays has been deter- 
mined so as to keep the watt consumption of the relays permanently on a level 
of 5,5 watts 


Battery voltage Resistance of the feeder | Max. loop re- 
bridge sistance 
24 volts 2x (40+ 25) 670 
36 volts 2x (90+ 60) 900 
48 volts 2 x (150 +-100) 1100 
60 volts 2 x (240 + 160) 1200 


With the increase of the number of telephone exchanges the average 
length of subscriber’s loops is gradually decreasing. Since for reasons of stress 
it is unlikely that diameters below 0,4 millimetre should be chosen for wires, 
a loop of 1100 ohms represents a distance of about 3,5 kilometres. This length 
appears to be satisfactory, inasmuch as the upper limit of 0,45 N for the atte- 
nuation specified by CCITT is reached sooner than the loop resistance of 1100 


ohms. 
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Consequently, as far as feeder bridges are concerned the 48-volt battery 
appears to be satisfactory and the 60-volt battery affords no particular advant- 
ages against the 48-volt type, not even on this count. 


2. Sparking 


The value of the voltage arising at opening an inductively loaded contact 
may be written as 


Assuming a relay to perform a definite task and an identical upper limit 
for the consumption at varying voltages, the current flowing in a 60-volt 
relay will have to be by 20 per cent below that flowing in a 48-volt relay. 
On the other hand, to ensure a satisfactory pull force, or an appropriate 
number of ampere turns, for a 60-volt battery the number of turns will have 
to be increased by the same percentage. 

While the value of di/dt will slightly decrease (however, not linearly, 
as the current disappearing period dt decreases almost linearly) with the cur- 
rent intensity, the self-induction will increase proportionally to the value 
of n*, if a 60-volt battery is used. At identical operating conditions, therefore, 
sparking will become more intense in case of a higher battery voltage. 

Note that in inductively loaded circuits the potential drop across the 
condensers will increase linearly with the battery voltage, e.g. when spark 
quenchers are used. Consequently, when condensers of the same make are 
used the safety factor, too, will drop by about twenty per cent, when the 
battery voltage is 60 volts. 


3. Battery and busbars 


As far as storage batteries and battery chargers are concerned the battery 
voltage has but little to say. The costs of the power supply plant represent 
a low percentage only of the overall costs of telephone equipment. Accordingly, 
slight fluctuations of prices, in one sense or the other, are for all practical 
purposes negligible. For 60 volts a larger number of cells may be required, 
while for 48 volts the dimensions of the cells may be larger. 

Busbars are in general dimensioned so as to make allowance for a certain 
voltage drop of let us say one per cent at the maximum load. Accordingly, 
provided the watt consumption remains identical, when a 60-volt battery is 
installed instead of the 48-volt, current intensity will drop by 20 per cent, 
and the cross-section of the copper bar by a further 20 per cent, making a total 
of 36 per cent. However, this reduction of weight will at the usual battery- 
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equipment distances amount to not more than 100 to 120 kilogrammes for two 
busbars and for a 10 000-line city exchange. Copper requirements for an ex- 
change of this size amount to about 20 000 kilogrammes approximately, one 
third whereof being accounted for by the windings. The amount of copper 
converted into busbars is too insignificant as compared with the total quantity 
of copper absorbed by the exchange. 


4. References 


The fact that the majority of crossbar exchanges operated all over the 
world uses 48-volt batteries, e.g. in the USA, cannot be ignored. If 60-volt. 
batteries gave any particular advantages, it may be taken for granted that 
60-volt batteries would have been adopted for the Panel System, instead of 
the 48-volt type, the more so as 60-volt batteries are employed in Strowger 
equipment from the earliest date. Notwithstanding this fact, the 48-volt battery 
continues to be in use in the various types of crossbar equipment. 

Ericsson crossbar exchanges are in general operated from 48-volt batte- 
ries, and only in exceptional cases designing engineers deferred to particular 
wishes expressed by certain administrations, e. g. in Finland, where the cross- 
bar exchange was designed for operation from 60-volt 5atteries. 


V. Conclusions 


The logical conclusion, from what has been set forth in this paper, is 
that in crossbar equipment 48-volt batteries have certain definite advantages 
as compared with 60-volt batteries. In the Strowger System the use of 60-volt 
batteries was necessitated by the step-by-step type operating magnets, rather 
than by any specific advantages of the 60 volts. 

It would be an exaggeration to state as if the problem of voltage were 
a determining factor in the cost estimates of an exchange. On the other hand, 
no conclusive reasons may be cited substantiating the use of higher voltages. 

The entire problem of voltage may have to be considered from altogether 
different aspects, when electronic elements will be introduced in the control 
system of crossbar exchanges. Then the necessity may arise of various power 
supply systems, and the magnets of the contact field may have to be made 
subject to the requirements of the controlling elements. However, as long as 
predominantly electromagnetic devices are used for controlling purposes, 
and other, e. g. not technical considerations, do not influence the decision of 
the designing engineer on the question of battery voltage, it is beyond doubt 
that 48-volt batteries will win the day against the 60-volt type. 
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Summary 


In development work of crossbar telephone exchanges the problem of the operating 
battery voltage for this type of equipment is considered. The paper gives a survey of some 
of the aspects which may havea bearing on the ultimate decision in the matter of voltage, 
and on hand of a number of points considered concludes in favour of the 48-volt battery. 


Professor L. Kozma, Polytechnic University, Budapest XI., Stoczek utca 2., 
Hungary. 
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UBER DIE NUMERISCHE LOSUNG DER MIT DEM LAPLACE- 
SCHEN AUSDRUCK ZUSAMMENHANGENDEN 
DIFFERENTIALGLEICHUNGEN 


Von 
G. Fopor 


Lehrstuhl fiir Theoretische Elektrotechnik an der Technischen Universitat, Budapest 


(Eingegangen am 5. September 1958) 


Bei der Behandlung technischer und physikalischer Probleme taucht 
haufig die Notwendigkeit der Lésung partieller Differentialgleichungen auf, 
die den Lapuaceschen Differentialoperator A enthalten. Es sind dies die 
Lapuacesche Differentialgleichung : 


AU = 0, (1a) 
die LapLtace—Porssonsche Differentialgleichung : 


AL =". (1b) 
die Wellengleichung : 
AU —k?U=0, (1c) 


und die inhomogene Wellengleichung : 
AU —K#U=S. (2) 


S bezeichnet in diesen Gleichungen eine gegebene Ortsfunktion, k ist im 
allgemeinen eine bekannte Konstante, kann aber auch eine bekannte Orts- 
funktion sein. U ist die gesuchte Funktion, die wir Potential nennen werden, 
obwohl sie in physikalischem Sinn oft kein Potential ist. 

Fir die Lésung der Gleichungen stehen uns in gewissen Sonderfallen 
mathematische Hilfsmittel zur Verfiigung, deren Anwendung jedoch meistens 
groBen mathematischen Aufwand erfordert. Die weitere Schwierigkeit besteht 
darin, daB die Lésung im allgemeinen eine unendliche Reihe oder ein Integral 
ergibt, deren Auswertung langwierig ist, und oft einen eigenen Apparat bean- 
sprucht. 

Der Gedanke liegt nahe, das Problem bei gegebener Geometrie numerisch 
zu lésen. Diese Methode hat den grundsatzlichen Nachteil, da8 die Anderung 
der Abmessungen dem Wesen nach eine neue Aufgabe stellt. Dementsprechend 
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ist die Ermittlung einer allgemeinen (fiir alle Abmessungen gultigen) Tabelle 
oder eines Schaubildes nur nach langwierigen numerischen Berechnungen 
méglich. Das Verfahren hat jedoch den grofen Vorteil, daB seine Anwendung 
keine groBen mathematischen Vorkenntnisse erfordert. Im konkreten Fall 
fiihrt die numerische Methode schneller zum Ziel als die analytische. In den 
meisten Fallen benétigt die Auswertung der allgemeinen analytischen Lésung 
fast so viel numerische Arbeit wie das unmittelbare numerische Verfahren. 
Daher bedeutet die »geschlossene Formel« oft ein triigerisches Ergebnis. 

Im weiteren sollen die Méglichkeiten der numerischen Lésung der 
erwahnten Differentialgleichungen in einigen Koordinatensystemen unter 
der Voraussetzung untersucht werden, daB die Funktion U an der Bereichs- 
grenze gegeben ist (erstes Randwertproblem). Wir spezialisieren die Aufgabe, 
indem wir nur Falle untersuchen, bei denen das Potential bloB von zwei 
Veranderlichen abhangig ist. In praktischen Fallen — infolge der Symmetrie 
der Anordnung oder anderer Ursachen — ist diese Bedingung recht haufig 
erfillt. Zuerst soll die bekannte Methode fiir die Lésung der ebenen LAPLACE- 
schen Differentialgleichung in rechtwinkligen Koordinaten (Verfahren von 
LIEBMANN) besprochen werden, worauf die Lésung der inhomogenen Wellen- 
gleichung (2) folgen soll, da die tibrigen Gleichungen als deren Spezialfalle 
auffaBbar sind. Der Gegenstand unserer Untersuchung ist die Ermittlung 
einer neuen, unabhangigen oder abhangigen Veranderlichen, die die numerische 
Arbeit erleichtert. 


1. Die Laplacesche Gleichung in der Ebene in rechtwinkligen Koordinaten 
[1, 2, 3, 4] 


Die bekannteste der numerischen Lésungsmethoden ist diejenige fiir 
die Lésung der ebenen Lapiaceschen Gleichung, das sogenannte LIEBMANN- 
sche Verfahren. Es besteht im Wesen aus folgendem : 

Es wird die Lisung der ebenen LAPpLAcEschen Gleichung 


0? U 27620 0 
0 x2 Aye S 


gesucht, wenn der Wert des Potentials U lings einer (geschlossenen) Kurve 
gegeben ist. Versehen wir den Bereich mit einem gleichmaBigen rechtwinkeli- 
gen Gitter, wobei der Gitterabstand d viel kleiner ist als irgendeine lineare 
Abmessung des Bereiches (Bild 1). Untersuchen wir den Zusammenhang zwi- 
schen dem Potential eines beliebigen mit 0 gekennzeichneten Punktes und 
den Potentialwerten der benachbarten Gitterpunkte (1, 2, 3, 4). Ist der Gitter- 
abstand d geniigend klein — wie dies soeben vorausgesetzt wurde —, lassen 
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sich die Anderungen mit einem Taytorschen Polinom dritten Grades gut 
annahern und es wird 


Addiert man diese vier Gleichungen, gelangt man zu der Beziehung : 


aU 


82U) d?  38U) d3 
om tgs (Ot # tu) a 
: x OX} Ox}, 2 eee 
ry ree P pas d? ay d3 
Ce oo wy ly 2 1 8053) 
en. aU) d@ (aU) @ 
; M G5 Oxy Gy 2 Salat 
eg a8 a eu BAP! a 
tc decelerate faye], 2 ae 


Ug Ue UU, Up 4D od? 


opal by 


8x? 


8 y? 


027U 


0 


(4) 


(5) 


Die Lariacesche Gleichung in der Ebene (3) ist erfiillt, wenn der Klammerwert 
gleich Null ist. Die Bedingung hierfir lautet : 


1 
EE etre cote U, + Us; + U4). 


(6) 


Die die Lapuacesche Gleichung befriedigende Gitterpunkt-Potentialverteilung 
ist somit gefunden, wenn der Potentialwert jedes Gitterpunktes das arithmetische 


Mittel des Potentialwertes der vier benachbarten Gitterpunkte darstellt. 


Bild. 1. Das Netz im rechtwinkligen Koordinatensystem 


Am Rande eines krummlinigen Bereiches sind die Gitterabstande d 


nicht samtlich gleich (siehe die rechte Halfte in Bild 1). In einem solchen Fall 


4* 
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schreiben sich die Anderungen, mit einem Taytorschen Polinom zweiten 
Grades ausgedriickt, zu 
362U 


8 x? 


2 
Bs 
? 


aU 
~U —| d,+ 
U; re eval : i 


fs) 
U,~ u+ (5 | d+ | 
9¥ Jo 


6*U| dz 
O47] 5 De . 
(7) 


aU 32U)\ d2 
Us Ue— [Fo] at al 
0 


aa2}, 2° 

a7 U. \ .d2 

U,~U, ee d, + | ee 
9¥ Jo Oy" Jo 2 


Diese vier Gleichungen ergeben zusammen mit (3) fiir die fiinf unbekannten 
GréBen im Punkt 0 (U, und die vier partiellen Derivierten) gerade finf lineare 
Gleichungen, deren Lésung lautet : 


0 


et dy d, ds d, | 1 = | 1 & ral 


= -|— ————- (8) 
d,d; + d,d,\d, + d3\d, ds d, +d, \d, dy 


Das Verfahren setzt sich sodann folgendermaBen fort : Am Rande des 
Bereiches ist das Potential unserer Bedingung gema® bekannt. Die Aufgabe 
besteht darin, das Potential der einzelnen Gitterpunkte derart zu bestimmen, 
daB diese die Gleichungen (6) bzw. (8) befriedigen. Praktisch geschieht dies, 
indem man die Potentialwerte — dem Gefihl nach richtig, sonst aber beliebig— 
annimmt und die so angenommenen Potentiale gemaB (6) baw. (8) korrigiert, 
wobei man alle Gitterpunkte systematisch durchgeht. Das Verfahren ist 
mehrmals und so lange zu wiederholen, bis sich innerhalb der Rechengenauig- 
keit keine weitere Korrektion ergibt. Sind die Potentialwerte der Gitterpunkte 
bekannt, kénnen die Aquipotentialkurven leicht konstruiert werden. 

Die numerische Arbeit 1a8t sich mit gewissen Kunstgriffen iibersicht- 
licher gestalten (»Methode der Reste«). Es mag erwaihnt werden, da8 man sich 
zur Beschleunigung der Konvergenz anfanglich zweckmaBig zu einer Uber- 
Korrektion der Potentialwerte entschlieSt. Es kann iibrigens nachgewiesen 
werden, dafs die urspriinglich angenommenen Potentialwerte die Konvergenz 
héchstens verzégern kénnen, da jedoch das Verfahren bestimmt zu richtigen 
Ergebnissen fiihrt. Das bedeutet zugleich, daB ein im Zuge der Berechnung 
eventuell vorkommender Rechenfehler automatisch korrigiert wird. 
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2. Die allgemeine Lésung des zweidimensionalen Problems 
Als allgemeinsten Fall soll nun die Differentialgleichung 
AU —k?U =S (9) 


untersucht werden, die die inhomogene Wellengleichung darstellt. Ist S = 0, 
gelangt man zur homogenen Wellengleichung, ist k = 0, zur Porssonschen 
Gleichung und, sofern k = 0 und S = 0, zur LariaceEschen Gleichung. 

Anderseits soll die Aufgabe derart eingeengt werden, daB die Unter- 
suchungen auf Falle beschrinkt bleiben, bei denen die GréBen nur von zwei 
Koordinaten abhingen, d.h. 


U =U (245%), S = S (x15 %9)- (10) 


Bedeuten x,, x, und x, allgemeine orthogonale Koordinaten, dann lautet 
bekanntlich die allgemeine Form des Lapiaceschen Ausdruckes [4] : 


AU = 1 le Belt 0 ae ) pace 
818283 |9%1! 81 Ox Ox, | 2 9X2 8x3 | 83 9%, 


| (11) 


Hier sind g,, g, und g, die MaBfunktionen, die durch den nachstehenden Aus- 
druck des Quadrates des elementaren Abstandes definiert sind : 


ds? = (g, dx;)? + (g2.d%2)? + (83 ds)”. 


Im gewahlten Sonderfall nimmt die Differentialgleichung (9) folgende Form an : 


1 & 8283 eae 0 (8381 alan (12) 
818283 1 9%1 | 81 0 xy Ox, | go 9x2 
Nach dem Differentieren wird 
126° 1 ) Le oU ot 
gi Oxi 818283 9% \ 81 0x 
1 0U 1 fe) leas (13) 
ge 0x2 81 8283 O%2\ Se ) 9% 
Fiihrt man die Abkiirzungen 
{ 
1 S) Ed — Gi i fs) ba = Gs (14) 
§18283 Ox1\ 81 818283 9%2 \ Se 


ein, nimmt die zu lésende Differentialgleichung folgende Form an: 


ple Ug OU wus: (15) 


e 0 x? Ox, 43 0 x2 0X, 
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Bild 2 zeigt einen Teil des zweidimensionalen Netzes. Die a-Werte — 
bedeuten Koordinaten-Differenzen und stellen nicht unbedingt einen Abstand 
dar (z. B. Winkel-Differenzen). Bei Auftragung des Netzes ist darauf zu achten, 
daB die Gitterabstinde im Verhdltnis zu den linearen Abmessungen des Berei- 
ches klein sein miissen, damit sich das Potential U (und die Funktion S) 
zwischen zwei benachbarten Gitterpunkten nicht sehr andere. Behandelt man 


Bild 2. Das Netzteil in allgemeinen orthogonalen Koordinaten 


eine Wellengleichung, so muB der Gitterabstand kleiner sein als der Parameter 
1/k, der eine Lange als Dimension hat, es mu8 also auch die Bedingung 
k g,a; <1 erfillt sein. 

Werden die Anderungen durch ein TaAyLorsches Polynom zweiten Gra- 
des angenahert, ergeben sich nachstehende Gleichungen : 


aU 2 2 
U, x U,+ |—} a, E te 

Xho ON lene 

9U (a2 2 
U, ~ U, 4 A> |: oy 

Xo} 9 O Nee 

(16) 

Uy Uy— | | a 

9X4) Ox 0 Pe 

aU 2 
U,~ U,— |_| ay cau oh 

X59 9 x25 2 


Gleichung (15) kann fiir den Punkt 0 in folgender Form aufgeschrieben werden: 


aU 


al 


oU 


A 


‘ 2 
0 810 


So= —k*U, + Gio t+ Goo 


0 


(17) 


salt 1 fae 


2 
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Betrachtet man U,, U,, Us, U, und S als bekannte Werte, lassen sich die finf 
Unbekannten ermitteln, doch benétigen wir nur U,: 


M 
U)= N ’ (18) 
M=U, 42%! : +% 6) +0, 9F [- ee Cao) + 
a, Sio 2 A (8% 2 
14 1 
+0, — 6,9] +0. =%| = Gy] —"P (as +) (ta + 0 
as Sio 2 oe 830 2 2 (19) 
Mae ce et : eee) Gro ee : +4 Gr] + 
Oy Sio 2 uo) §20 2 
A> | a a Gard a k2 
25 as ‘| pe ‘6,4 + ae Peat 2 Ga) +5 (on + 25) (ta + 8) 
Os Si0 2 “a 820 : 2 


Sind die denselben Koordinaten zugehérigen Gitterabstande einander gleich, 
wird der Ausdruck etwas einfacher; 4, = 43 sowie Gg = dy gesetzt, gilt 


a 1 a af al a. a 1 a 
M =U, (+2 Cy] + Use 4% Gy] + Uo (—- — Cra 4 
& \ S10 2 A, \ $30 2 a \ S10 2 
a,{1 a 
+U,+ = Ge — So % o>» (20) 
a, \ $20 2 
N==2° : pes eA crore. 
A, Zio A, 830 


Bei Untersuchung der verschiedenen Koordinatensysteme kann man 
nun derart verfahren, da®8 man mit den bekannten MaBfunktionen g,, 82 83 
die Beziehungen G, und G, gemaB (14) bildet ; setzt man diese in Gleichung (19) 
bzw. (20) ein, stehen die Zusammenhange fiir die Korrektion des Potentials 


zur Vertiigung. 


3. Kartesische Koordinaten 
[1, 2, 3, 4] 


In rechtwinkligen kartesischen Koordinaten gilt 


%==—s=—1, (21) 
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woraus G, = Gz = 0. (19) baw. (20) nehmen mit a; = d; die Form 


d,+d, S 
py eth Us, d, +d; ae dy +dq utto = 3 + dy) (dy + dy) 
dy d, d. 4 


dtd at+ds td dg th Mg aya, +d, 
Ge. d, Gs a $ 


= 
(22) 


an (hier bedeuten die Gitterabstinde wirkliche Abstande), bzw. falls d, = 
— dy = d, = d, = d, wird 


_U,+U, + U, + Uy — Sod? 
4 + k2d? 


U, (23) 


Ist S = 0 und k = 0 dann geht (22) in die unmittelbar abgeleitete Gleichung 
(8), bzw. (23) in die Gleichung (6) iiber, wie dies auch sein muB. 


4. Polare Koordinaten in der Ebene 


Aus den Ausfihrungen unter 1. und 3. ist ersichtlich, daB sich fiir einen 
im Inneren des Bereiches liegenden Punkt ein einfacher Zusammenhang fiir die 
Potentialkorrektion ergibt. Dagegen ist die Ermittlung des gewogenen Mittel- 
wertes fiir einen am Rande des Bereiches liegenden Punkt sehr langwierig. 
Da in der Praxis der kreisférmige Umrif haufig vorkommt, lohnt es sich, 
die Lésungsmethode auch fiir polare Koordinaten zu untersuchen. 


a) Korrektionsbeziehung in den urspriinglichen Koordinaten [3] 


Das polare Koordinatensystem ist nichts anderes als ein Zylinderkoor- 
dinatensystem, in welchem das Feld unabhingig von der axialen Richtung 
betrachtet werden kann (Problem in der Ebene). Die Koordinaten und die 
MaBfunktionen sind 


y= T, X= Py, Xe=—Bz 

f= 1, B= 5, g=1; (24) 
1 

G6, =—— (r) =—, G,=0 
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Bezeichnet man die radiale Unterteilung mit d, die azimutale dagegen 


mit a (Bild 3), folgt aus (19), daB 


Bild 3. Das Gittersystem in polaren Koordinaten 


plea (eee 14 4 ut AS + 
1 21 a, 73 d, 27 
d S 
4 UBS _ 29 (a, + ds) (ta +)» (25) 
a,re 2 
2 
pelea is he d, hte Foe BS) 4S + dle t a 


1 


Ist-d, =d,=—d und 4,—4,=4, 80 gilt gemaB (20) 


a*To ¢ (2 6) 


- (27) 
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einzufiihren, der bei gleichmaBiger radialer Teilung die von innen gerechnete 
Ordnungszahl des betreffenden Kreises bedeutet. Damit wird 


1 
p+ s—]ui+ [1-2] tet age + Ud Som 
in 2n 2n n?a (28) 
1 
2 [I+ eS 4+ kd? 
n*~a 


Gesondert muB der Fall untersucht werden, wenn auch der Punkt r = 0 
in dem betrachteten Bereich liegt. Der Mittelpunkt hat selbstverstandlich 
viele benachbarte Punkte, und infolge der Singularitaét der Differential- 
gleichung gelten dort auch die Korrektionszusammenhange nicht. Da samt- 
liche Gitterpunkte vom Mittelpunkt gleich weit entfernt liegen, ertibrigt sich 
jetzt ein Abwagen. Legt man vom Mittelpunkt in verschiedene Richtungen ein 
rechtwinkliges Netz, lat sich (23) folgendermaBen verallgemeinern : 


CS :Tr=0, n=0; (29) 


wobei U; das arithmetische Mittel des Potentials des m benachbarten Gitter- 
punktes bezeichnet : 


Ups ee SU (30) 


b) Einfiihrung einer neuen unabhdngigen Verdnderlichen [3] 


Wie man sieht, gestaltet sich die Berechnung dieser Zusammenhange 
viel umstiandlicher als derjenigen im kartesischen Koordinatensystem. Die 
Schwierigkeit ergibt sich dadurch, daB die Funktion G, von Null verschieden 
ist. Man kann versuchen, dies so auszuschalten, da8 man den radialen Gitter- 
abstand nicht als gleichbleibend annimmt, was im wesentlichen so viel bedeu- 


tet, daB man statt der Verdnderlichen r eine geeignete, neue unabhangige 
Verinderliche 0 (r) einfiihrt. 


Die neuen Verdnderlichen sind 


$1 9” S2=T(0), g,=1; (31) 


ea cael 
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Aus der Bedingung G, = 0 folgt, da® die in der Klammer stehende Funktion 
konstant sein muB. Der Wert der Konstanten wird der Einfachkeit halber 


gleich 1 gewahlt, womit sich die fiir die Ermittlung der Funktion r (g) dienende 
Differentialgleichung in der Form 


i 
T = r(e) (32) 
schreibt. Das Integral dieser Differentialgleichung lautet 
r 
r=Re®, o=In—, 33 
Q cs (33) 


wobei R eine unbestimmte Konstante mit einer Linge als Dimension ist. 


Die Ma8funktion g, hat den Wert 
gp =— =r=Re. (34) 


Bezeichnet man die Gitterabstande in der Richtung y mit a, in der Richtung 
o dagegen mit f, gilt mit G, = 0 aus Gleichung (20) 


U, fe +0, + Uae Po ais 26 
oe To aro a TS (35) 
a 
pA - 2 ke 
pois eee af 


Zur einfachsten Beziehung gelangt man, wenn man a = f wahlt ; in diesem 


Fall gilt in weitgehender Analogie mit (23) 
= Uy U, + Ut Uy paren 


U 36) 
: 4+ hk? ra? ( 
Sind die Gitterabstinde ungleich, wird gemaf (19) 
Yaa ot Sain | § eet! 13 yy, pes au “4U, (Rie fama 5 
oF ar 3 Oy 
50 72 37 
Se (2 aa) (a 0) (37) 


yoMte ate ate ety A 75 (ay + 45) (4g + %)- 


ay Ay Ag O4 


Die azimutale Unterteilung a wird auf Grund der Geometrie angenom- 
men. Durch zweckmaBige Wahl des Parameters R kann man f = a erreichen. 
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Bezeichnet man den gréBten Halbmesser mit rj, lassen sich die Halbmesser 
des Netzes der Reihenfolge nach aus der Beziehung 


= T)+1 es 
r,=Re®, rt, =Ree*; —H=e*%=e (38) 
T, 
ermitteln, and zwar 
9 =e A-l1 
Fie Ty, Ta = C815 hy = C-Fiot-2 ey 1 eee (39) 


Ein solches Netz zeigt Bild 4. Man erkennt, daB der radiale Gitterabstand 
nach auBen zu anwichst, was gleichzeitig auch den grundlegenden Nachteil 


Bild 4. Das Gittersystem in polaren Koordinaten mit den neuen unabhangigen Veranderlichen 
oe =In R/r. In diesem Fall ist a = 2/12 = 15° und daher c = e—¢ = 0,77 


dieser Methode beleuchtet : In Nahe der Achse haufen sich die Gitterpunkte 
sehr an, was die Rechenarbeit tiberfliissig anwachsen la8t ; am auBeren Rande 
des Bereiches verschlechtert sich gleichzeitig die Genauigkeit des Verfahrens. 


Dies verursacht besonders dann grofe Schwierigkeiten, wenn der Bereich 
auch den Punkt r = 0 enthalt. 


c) Einfiihrung einer neuen abhingigen Verénderlichen 


Ein dritter Weg zur Lésung kann folgender sein : Man fuhrt statt der 


Veranderlichen U(r, @) eine neue abhingige Verdnderliche mit Hilfe der 
Gleichung 


U (r,9) =r v(r, @) (40) 
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ein. Der Exponent m hat vorerst einen beliebigen Wert. In der Differential- 
gleichung (12) wird sich der Wert des ersten Gliedes andern und es wird 


1 fs) ee aU) _ oe aaa le 
81828;9%1\ 81 9% r Or Or 
(41) 
2 
as phe + (2m + 1) pea oy mtr" vy, 
is or 
1 
Ist m= — ie verschwindet der Faktor des ersten Derivierten, was formal 
bedeutet, als ware G, gleich Null. (13) nimmt die Form 
on 2 1 02 n 
ar + m2r™ py + —r™— — Br™v=S; 42 
Orns F250 gi? a 
an, wobei 
1 
ist. Vereinfacht man durch r™ und substituiert man m = — >; So gilt 
oe 2 - 
sib Saye ome _ (x2 |» = 8 ]r- (44) 
Oye tr ne 4 Ar? 
Nach Vergleich mit Gleichung (15) kann man formal schreiben, da8B 
——-> 1, G6, 95 g2-—>- I> G, —— 0, 
&1 1 2 (45) 
k2 > ke? : tS St 
4 r? 


wobei r an der Stelle zu nehmen ist, wo die Korrektion geschieht, d. h. r = 7. 
Bezeichnet man den radialen Gitterabstand mit d, den azimutalen 
dagegen mit a, wird aus (20) im Falle gleicher Gitterabstande 


d? — 1, 
V1 + Us Sere (vg + V4) — So Vr d° 
a 


= ss 46 
Vo 9 1 re d2 | i 2 re aa op ( ) 
| Ors | Are 
Fiihrt man wieder die Ordnungszahl der Kreise n = r,/d ein, gilt 
1 a 
Wires chee g (vq + V4) — So/n Vas 
n* a 

v= (47) 


1 1 


n2a2 8n? 


4 ka? 


2 (1+ 


« 
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Allgemeine Gitterabstinde vorausgesetzt, ergibt sich aus (19) 


d,-+d Sy y— 
M= be Ba = “8b, Sie web j= met At, peat Say, ak ee Vr (d,+d3) (a2+44)+ 


2 
1 A, 3 A479 


na tettey tds te dod dir ds 


2 2 2 
d, Gs Fé d, a,r2 8 r2 


(42+ a) + (48) 


k2 
ci = (d, + dg) (42 + a4). 


Wieder mu8 der Fall separat untersucht werden, daB das Potential 
auch im Punkte r = 0 ermittelt werden soll. Da das Potential U und seine 
Derivierten iiberall endlich sind, lassen sich folgende Beziehungen schreiben : 


v= )rU—0 falls r—0, 
Oo ee a U— us U=——, falls Oe (49) 
Or oir 2\r 2\r 2r 


Im Punkt r = 0 streben also simtliche Derivierte von v auf r gegen unendlich. 
Dem Halbmesser r = d (n = 1) entlang kann also nicht v, = 0 substituiert 
werden, da dort die Annaherung mit dem Taytorschen Polynom zweiten 
Grades nicht gestattet ist. Die fiir die Ermitthung der Korrektionsbeziehung 
nétige vierte Angabe ist eben durch (49) gegeben. An der Stelle r = d lassen 
sich namlich folgende Beziehungen aufschreiben : 


d? 
/0 2 


82» 


Or 


ney t |: | d+ 
0 


r 


vo 
7 Vg + ——d 
Bea 5 


(50) 
2 v 


8 gy? 


a2 


° 


0 2 

2 2 
«+f ar 

0 Oy), 2 


Bericksichtigt man nun auch die auf den Punkt r = d bezogene Differential- 
gleichung (44), gelangt man zu der Korrektionsbeziehung 


nme + || a+ | 
0 


Ov 
Oy 


nwt | 


1 ae 
20+ ae ee Ps) — Sy \/a? 


OS 9 , neal 50 
2,75 + — + k2 d2 
a2 
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Vergleicht man die Ausdriicke (28) bzw. (47), 1a8t sich der Vorteil dieser 
Methode feststellen : Wahrend namlich die Potentiale U, und U3; gewogen 
werden miissen, ist dies bei den Hilfspotentialen v, und vg nicht notwendig. 

Man verfahrt wie folgt. Vorerst wird das Netz angefertigt. Anstatt der 
- Randwerte schreibt man die Werte von v = \r U ein. Danach wird die 
~ Potentialkorrektur an Hand der Gleichungen (47), (51) bzw. (48) durch- 
gefiihrt. (Ist n > 3, bedeutet die Vernachlassigung des Gliedes 1/8 n? einen 
Fehler, der kleiner ist als 1%.) Aus den so erhaltenen Werten des Hilfspoten- 
tials v kann man mit Hilfe der Beziehung U = v/|/r die Werte des Potentials 

U errechnen. Das Potential des Mittelpunktes 1aBt sich aus (29) ermitteln. 


5. Zylindersymmetrische Felder in Zylinderkoordinaten 


a) Korrektionsbeziehung in den urspriinglichen Koordinaten (2) 

In zylindersymmetrischen Feldern hangt das Potential U nicht vom 
Winkel y ab. Die Koordinaten und die MaBfunktionen werden zweckmafig 
in folgender Weise geordnet 


gil, Se=—1., g=F3 (52) 
ik =} il 

Gia= 0 C= SST SS 
r or r 


Bezeichnet man die Gitterabstande in der Richtung mit d,, dz, in der Rich- 
tung r mit d,, d,, ergibt sich aus (19) 


d,+dy(,, 4 
ds 


27 


d 
M—U,H=" +0, ee Nees each 
2 


3 


pu PEA — 4\— Sat alae +40, 


2r, 
4 0 (53) 
vahth SES Bh: _ 4 Rte Fah _ 4, | ne 
d, d, 2 To ds d, 2 To 
2 
Bet (d, + ds) (dz + da) - 
Im Inneren des Bereiches gilt im Falle gleicher Gitterabstande aus (20) 
1 1 
U,+ (1+ |U.+ 0s + — 5 ,|Us— Soe 
Ue = sg ; (54) 


44+ kd? 
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wobei n = r,/d. Langs der Achse — wo r = 0 — gilt aus Symmetriegriinden 
U, = U,, woraus 
3 4+ kd? 


b) Einfiihrung einer neuen unabhdngigen Verdnderlichen 


Die Anwendung der Beziehung (54) gestaltet sich wieder deshalb schwie- 
rig, weil die Potentialwerte im Zahler gewogen sein missen. Das System (52) 
weist eine weitgehende Analogie mit dem System (24) auf, man kann also ver- 
suchen, die neue unabhangige Veranderliche 9 = Inr/R einzufiihren. Ganz 
dbnlich, wie in Punkt 4 wird 


r 
Ls 5 Rags Ore ee Ky =P 


Bu No 8a Rel nea eee (56) 
Peo e he 
G,=0, G, = ee Sey | 
Re? 3 0 Re? 


Bezeichnet d den Gitterabstand in Richtung z und f in Richtung a, ergibt 
sich aus (20) 
a2 


Oye Us A 
Op ro 8 


2 (1+ 


; (Us + Uy) — Spd? 


(57) 


ge 
# awe 

rp) 

(57) ist nicht einfacher als die Beziehung (54). Da ferner auch eine Ungenauig- 


keit wegen der ungleichen Gitterabstinde auftritt, ist es hier unzweckmaBig, 
diese Methode anzuwenden. 


c) Hinfiihrung einer neuen abhingigen Verdnderlichen 


Ks kann das Verfahren der Einfiihrung einer neuen abhangigen Verander- 
lichen angewendet werden. Da die zu Null zu machende Funktion G, aus dem- 
selben Glied stammt wie im vorigen Punkt, kann man die Ergebnisse von 
dort tibernehmen (Gleichungen 40—42) : 


ti(eae es v (z, r) (58) 


it, Se (59) 
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Bei gleichen Gitterabstanden folgt aus (20) 


gy ote (60) 
Ane aoe PES Ie 
4n? e 


Sind dagegen die Gitterabstanden allgemein, dann wird gemaB (19) 


M = Ste =e ott 4 gee nets 
1 2 3 4 
ge 
a itie (d, + dg) (dg + dy), (61) 


N 


_ ay + dy d, + dg d, +d, ds+d, (dy +s) (dy + 44) 
oe es oe, 8 r2 35 


k2 
oxo (d, + ds) (dz + d,)- 


Auch hier mu8 wieder der Fall gesondert behandelt werden, daB das 
Potential auch langs der Symmetrieachse (r = 0) zu bestimmen ist. Wie im 
vorigen Punkt (Gleichung 49) 1aBt sich auch hier zeigen, daB die partiellen 
Derivierten der Funktion v auf r im Punkt r = 0 gegen unendlich streben. 
Mit derselben Methode laBt sich auch die Korrektionsbeziehung an der ersten 
Teilung ermitteln. Schreibt man in (51) statt a formal 1 und paSt man die 
Indizes den vorherigen an, dann ergibt sich 


_ M+ 2 0, + 05 — So Vai n=l. (62) 
4,75 + k? d? 


Das Verfahren und die Vorteile der Methode sind die gleichen wie bei 
den Polarkoordinaten besprochen, ja hier brauchen sogar die Werte des Hilfs- 
potentials nicht erwogen zu werden. Verzichtet man auf eine gréBere Genauig- 
keit als 1%, kann in (60) bei n > 2 das Glied 1/4 n? neben 4 vernachlassigt 
werden. Ist k = 0 (Lapiacesche Gleichung, Porssonsche Gleichung), erleich- 
tert man sich die Rechenarbeit sehr mittels einer Tabelle, die folgende Spalten 
enthialt 


Vo 


n 1 2 >2 
x x x 


4,75 3,9375 4 


Die Tabellenwerte lassen sich mit der Rechenmaschine schnell ermitteln. 
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6. Zylindersymmetrische Felder in Kugelkoordinaten 


a) Korrektionsbeziehung in den urspriinglichen Koordinaten 


Wegen der Form des Umrisses ist es haufig zweckmabBiger, die zylinder- — 


symmetrischen Felder in Kugelkoordinaten zu berechnen. Infolge der Zylinder- 
symmetrie ist das Potential U unabhangig vom Winkel ~. Dementsprechend 
ergeben sich folgende Koordinaten und Maffunktionen : 


Coif ae Kan te = Gs 


(gyotk- i—i- g3=rsind; 


= I oe aan ee 
r2sin? Or r 
(63) 
ne 1 nog <j __ cos 0 
>" 2sin 8 OB r2sin 3 


Bezeichnet d den radialen Gitterabstand und a denjenigen in Richtung des 
Winkels, dann wird mit dem Parameter n = 1,/d aus (20) 


m= (1+ —|t+35 [i+ Su + (1+ 
n n? a? 2 sin By n 
1 
sarees deere (64) 
na? 2 sin Jy 


N=2(l+— 3) + ee. 
| noe 


Da dieser Zusammenhang au8erordentlich umstandlich ist, kann er praktisch 
kaum angewendet werden. Sind die Gitterabstande ungleich, so ergibt sich 
sebstverstandlich eine noch uniibersichtlichere Formel. 


b) Einfiihrung einer neuen abhdngigen Verénderlichen 


Es kénnte versucht werden, die Vereinfachung mit Hilfe einer neuen 
unabhangigen Verdnderlichen zu erreichen, doch bedeutete dies keine groBe 
Vereinfachung, und die ungleichmaBigen Gitterabstande verursachten weitere 
Nachteile. Es soll daher eine neue, zweckm&Bige unabhangige Verdnderliche 
ay werden, bei der die Funktionen G, und G, einfacher, woméglich Null 
werden. 


Man fihrt die durch die Gleichung 


U (r, 3) = r“ sin’ 9 v (r, 0) (65) 
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definierte neue, abhangige Veranderliche v ein. Der Lapiacesche Ausdruck 
(11) nimmt dann folgende Form an: 


1 1 
U = ie preter te sin" Oe 42 |sin a2 psindv . (66) 
r?sin@ Or Or Od Od 
Fihrt man die Differentierungen durch, wird : 
resin’?! #[ , 020 F) 
AU = 7) 2 4 ou + rt w(t Dot 
r?sin ? Or? Or 
{ (67) 
6? v cosdOv . vcos*# — sin? d 
— 2y+1 vi. 
tag? | . in 080 sin? } | 
Die Differentialgleichung erhalt schlieBlich die Endform 
2 2 
au og ttle ay) ees LeeP Oa oe OE OUeey 
Or? re OF Me roy r2 sin? 0d 
' (68) 
veos?F—sin?D aa S 
r2 sin? ? r“ sin’ & 


Vergleicht man diesen Ausdruck mit der allgemeinen Form (15), so ergeben 
sich folgende Kongruenzen 


1 2 1 Oi) 
Gems Lg by Apri yy i ea ene as war 
r r2 sin? 
(69) 
Pe eee: 
eyo Shae Plee B M sin Ealc rai we > ¥ 
r r? sin? r¥ sin” 0 
Durch Wahl von » = — 1, bzw. y = — 1/2 lassen sich die Funktionen G, 


bzw. G, gleich Null setzen, so da man nach Einfihrung der durch die Glei- 
chung 


ro) (70) 


definierten neuen Veranderlichen v mit folgenden Parameterwerten zu rechnen 


hat: 


g,—— 1, G,—— 9, g,——T» G,— 0, S ——>SrJsin@. 


(71) 


os cos? } + sin? ? 


a3 ed 1h 
k2 ——-> k? + —— k? + : 
r2sin? ? 2r2%sin 3 


1 + sin? ? 
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Mit n = r,/d erhalt man aus (20) 


Pica ar Praia y (v2 + 4) — Som sin 8d* 
v= (72) 
1+ sin? Dy 4B 
4n? sin? J, 
Sind die Gitterabstande ungleich, folgt aus Gleichung (19) 
d,+d 
Mis p SES pe ae ee ee 
1 Ay 75 ds ee 
> =° rosin By (d, + ds) (4, + 4), 
(73) 
ae (acts Oa ese wes teal joie - 
d, Gy Te. d; a,rz 
1 + sin a 
ite fea nai a + ds) (43 + 44) ls is. (ds + ds) (ag + 44) - 


Die Frage der Geraden sin # = 0 (Symmetrieachse) und der folgenden 
Netzgeraden (# = a bzw. # = x—a) muB wieder separat betrachtet werden. 
Dav=r Vsin 3 U, streben die Derivierten von v nach # an der Stelle sin 8=0 


Bild 5. Das Gittersystem in Zylinderkoordinaten 


(d.h. # = 0 und # = a, wo v = 0) gegen Unendlich, genau so, wie im Falle 
der Polar- oder Zylinderkoordinaten. Hier kann ebenfalls die Anniherung 


a2 = rein bo aA! ow rU ee (74) 
ao” 2Vsnd 2Vsnd 2sind 


ae 
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benutzt werden. In der Umgebung eines an der Stelle 0 = a liegenden Punk- 


tes gilt 
9 2 2 
04% + Bal d+ oe! a 
Or Jo Or} 2 
9 2 2 
ram 0— |e dl ee Cah 
OT's Or), 2 
715 
5 02 v) a? 2) 
Vg Ug + == a —| — Fre, 
% +——— 2 Sites 
oT sind, (603), °2- 


Beriicksichtigt man auch die Gleichungen (68) und (71), ergibt sich an der 
Stelle #®—=a folgende Korrektionsbeziehung 


y+ % + — 5 ta Som \sin ad? 
Vg = 9 7 a 1 Behl 
24 + sin*a ee 


n? a2 2n? sin? a n’asin a 


“(Tey 
7— a 


Ist a < 0,25 ~ 14°, so ist mit einem Fehler unter 1% sin a ~ a und es wird 


vy + vg + 2 y, — Son sin ad® 
n?a a 
24 2 4+—-+ kd 
vi n?a a an?! 


Sind an denStellen # = a bzw. 0 = x—a die Werte des Hilfspotentials 
v und damit auch diejenigen des Potentials U bekannt, 14Bt sich das Potential 
lings der Symmetrieachse gem4B (64) ermitteln, Da hier aus Symmetriegriinden 


U, = U,, gilt 


fb + fp ———| 0+ aU, 5,0 . 
Uy= ~— i = 0 eee! 
2 } + ce + hee = 
n? a? 


Eine weitere Singularitat weist jener Mittelpunkt r = 0 auf, der viele 
benachbarte Gitterpunkte besitzt. Da jetzt v =r sind U, ist der Wert des 
Hilfspotentials im Punkte r= 0 gleich Null, und seine Ableitungen auf r 
bleiben endlich. An den Stellen r= d (d.h. n= 1) ist also Gleichung (72) 
anwendbar, und nur der Punkt r = 0 ist gesondert zu untersuchen. Will man 
die Potentiale samtlicher benachbarter Punkte beriicksichtigen, so miBte 
man ein ziemlich langwieriges Wiegen durchfiihren. Da am innersten Gitter- 


* 
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kreis die Anderung nach @ ziemlich klein ist, genugt es, diejenigen Potentiale 
zu beriicksichtigen, die den Werten 9 = 0, n/2 und x zugehéren. Stellt man 
sich im Mittelpunkt ein raumliches kartesisches Koordinatensystem vor, 
kann man auf Grund des Gesagten leicht die Beziehung 


frees 
6 


Us re ae (79) 
1+— hd? 


ableiten, wobei U; den arithmetischen Mittelwert der Potentiale der sechs 


benachbarten Gitterpunkte bedeutet. Im zylindersymmetrischen Fall stimmt 
das Potential des Punktes ? = = mit dem Potential der auf dem »Aquator« 


liegenden itibrigen drei Punkte iiberein, es wird somit 


ee 1 It 
= 9=—0 AU |0 = — 
U, ALA )+ | ; 


+U@=2)), (80) 


Wie ersichtlich, ist Gleichung (72) von viel einfacherer Struktur als 
Gleichung (64), dieses Verfahren erleichtert mithin tatsichlich die Rechen- 
arbeit. 


Zusammenfassung 
Es wurden die Methoden der numerischen Lésung der in der Form 
AU—FU=S 


geschriebenen partialen Differentialgleichungen fiir den Fall untersucht, daB die Funktionen 
nur von zwei Verdnderlichen abhingen, und das Potential der Randwerte gegeben ist. Es 
wurde festgestellt, daB in kartesichen, in zylindrischen und in Kugelkoordinaten die Rechen- 
arbeit durch die Einfiihrung einer geeigneten neuen abhangigen Veranderlichen vereinfacht 
werden kann. Das Verfahren hat immerhin den Nachteil, daB in singularen Punkten (wie 
z. B. im Mittelpunkt} und in den mit ihnen benachbarten Punkten andere Beziehungen fii 
die Potentialkorrektion gelten als in den Punkten allgemeiner Lage. Die Konvergenz, die 
Fehlerschranken usw. der Methode wurden nicht erértert, weil die Struktur der Differential- 
gleichungen fiir die neuen Veranderlichen dieselbe ist wie diejenige fiir die schon friiher ein- 
gehende Untersuchungen durchfiihrt worden sind. 

Die verschiedenen Methoden sind in der beiliegenden Tabelle .zusammengefaBt. Aus 
der Tabelle ist ebenfalls ersichtlich, daB die hier vorgeschlagene Methode der Einfiihrung einer 


neuen abhangigen Verinderlichen (oder eines Hilfspotentials) im allgemeinen zu einfacheren 
Korrektionsbeziehungen fiihrt. 
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; Bezeichnungen 
 D=tgo Verlustfaktor Uc Amplitude der Briickenaus- 
ig (t) Momentanwert des Anoden- gangsspannungs-Komponente, 
stromes die mit der Speisespannung in 
IA Anodengleichstrom Phase ist 
Ala Unterschied der Anodengleich- Us; Amplitude der Briickenaus- 
strome gangsspannungs-Komponente, 
Alan Unterschied der Anodengleich- deren Phase zur Speisespan- 
stréme, hervorgerufen durch die nung um 90° verschoben ist 
n-te Harmonische Un Amplitude der n-ten Harmoni- 
Alan Unterschied der Anodengleich- schen in der Briickenausgangs- 
stréme, hervorgerufen durch spannung 
samtliche Harmonische Uo Amplitude der Briickenspeise- 
ij —) —1 imaginiére Einheit spannung 
kn Relativwert der n-ten Harmo- uy, (t) Momentanwert der Wechsel- 
nischen spannung am 1. Gitter 
Leff effektive Induktivitat Us (t) Momentanwert der Wechsel- 
In Normalinduktivitat spannung am 3. Gitter 
Lx unbekannte Induktivitat U4 Sperrspannung des 1. Gitters 
n Ordnungszahl der Harmoni- Uso Sperrspannung des 3. Gitters 
schen ZN Normalimpedanz 
N=|Zy| Absolutwert der Normalimpe- zu messende Impedanz 
danz X = |Zx| Absolutwert der zu messenden 
S(t) Steilheit des Eingangsgitters Impedanz 
der Mischroéhre A Relative Impedanzabweichung 
Sc Mischsteilheit An MeBfehler der relativen Impe- 
Smax Maximalwert der Eingangsgit- danzabweichung 
ter-Steilheit Dn Phasenwinkel der Normalimpe- 
Sn n-te harmonische Komponente danz 
von S(t) Dx Phasenwinkel der zu messenden 
U Amplitude der Briickenaus- Impedanz 
gangsspannung gy = &x—®yn Phasenwinkelunterschied der zu 
vergleichenden Impedanzen 
Wo Eigenfrequenz der Spule 
Q relative Frequenz 


1. Einleitung 


Bei der Lisung verschiedener technischer und wissenschaftlicher Pro- 
bleme erweist es sich haufig als nétig, die verhdltnism4Big kleine relative Ande- 


rung irgendeiner Impedan 
die GréBe, sondern auch di 


z zu messen. In den meisten Fallen mu8 nicht nur 
e Richtung der Anderung festgestellt werden. Man 


benutzt zu diesem Zwecke am besten eine Wechselstrom-MeBbriicke mit 
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phasenempfindlichem Indikator, der dengroBen Vorzug hat, gegen Stérspan- 
nungen wie z. B. Netzbrummspannungen, die an der zu messenden Impedanz 
auftreten, sowie gegen die Harmonischen des Generators, der die Briicke 
speist, unempfindlich zu sein. Die allgemein gebrauchlichen Schaltungen 
enthalten einen Réhrengenerator mit einer oder mit zwei Réhren, einen zwei- 
stufigen Verstérker und einen phasenempfindlichen Detektor mit Doppel- 
diode. 

Im folgenden soll eine schaltungstechnisch neue MeBbriicke mit phasen- 
empfindlichem Indikator beschrieben werden, deren Lésung sehr einfach ist, 
und die nur zwei Réhren mit je zwei Steuergittern (Hexode, Heptode, Penta- 
grid usw.) bendtigt. Die Stabilitat der Schaltung wird durch symmetrischen 
Aufbau und negative Rickkopplung gesichert. 


2. Die Differentialbriicke [1], [2] 


Zwei Impedanzen kénnen mit einer Differenzialbriicke verglichen 
werden. 
Die zu messende Impedanz betragt : 


die Normalimpedanz hingegen 
Zy — N él PN . 


Die Ausgangsspannung der Differentialbriicke, die durch einen Generator 
mit der Spannung U, cos wt gespeist wird, betrigt : 


1,,Zy—Z 
“= = Neoswt. 
Zx+ Zn 


Die Ausgangsspannung kann in zwei Komponenten zerlegt werden, deren eine 
die parallele, mit U, in Phase, die andere hingegen — die senkrechte — 
gegen U, um 90° in Phase verschoben ist : 


u = U cos (wt — ~) = U cos at cos y + U sin wt sin py = U, cos wt + U,sin ot. 


Die parallele Komponente der Briickenausgangsspannung betragt 


2. 8 7 eae eae 
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die senkrechte Komponente dagegen 


xX 
2 (“si 
Lyla 


U,= Im —U, X_N =U, 
xt+4n | aaeiye 


= cos 
N Ne 


wobei der Phasenwinkelunterschied zwischen den beiden Impedanzen 
YP — @ 5 hoe @x . 


Weichen die Absolutwerte der beiden Impedanzen wenig voneinander ab, 
betragt also der Absolutwert der zu messenden Impedanz 


X=N(1+A4) 
wobei A-z1, 
so gilt mit guter Annaherung 
U.= : U; a 
2 1+ cosy 

und 

Ue = Js Uy Sel a . 

2 1+ cosp 


Gelangen auf kapazitivem oder auf induktivem Wege Stérspannungen 
vom Leitungsnetz an das MeSobjekt — was in der Praxis oft vorkommt —, 
werden diese vom Indikator mitgemessen. Einen weiteren Fehler verursacht 


Bild 1. Die Differentialbriicke 


der Umstand, daB die Briicke, die auch Reaktanzglieder enthalt, auf die vom 
Speisegenerator erzeugten Harmonischen nicht abgeglichen ist. Wird ein 
phasenempfindlicher Indikator benutzt, so gibt die Messung nicht nur die 
Richtung des Unterschiedes der beiden Impedanzen an, sie ist auch gegeniiber 
Stérspannungen weitgehend unempfindlich. 
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In der Schaltung gem4B Bild 2 wird die MeBbriicke von einem Oszillator 
gespeist. Die Ausgangsspannung der Briicke wird erst verstarkt und dann 
einem Phasendetektor [3] zugefiihrt, der seine Vergleichsspannung von dem 
die Briicke speisenden Oszillator erhalt. Die Verstarkerstufe zwischen Briicken- 
ausgang und Phasendetektor ist erforderlich, weil der Phasendetektor einen 


Bild 2. MeBbriicke mit phasenempfindlichem Indikator 


Phasen- 
detector 


kleinen Eingangswiderstand und eine geringe Spannungsempfindlichkeit 
besitzt. 


3. Ausgangsspannung der Differentialbriicke, wenn die Speisespannung 
Harmonische enthalt 


Wie bereits festgestellt, ist die Differentialbriicke abgeglichen, wenn 
sowohl die Absolutwerte, als auch die Phasenwinkel der zu messenden Impe- 
danzen gleich sind. Haben die zu messenden Impedanzen gleichen Charakter 
(z. B. Kondensatoren mit in Reihe geschalteten Verlustwiderstanden), kénnen 
die Abgleichbedingungen bei jeder beliebigen Frequenz erfillt werden. Befin- 
det sich jedoch in dem einen Zweig der Briicke ein Kondensator mit in Reihe 


Bild 3. Messung verlustbehafteter Kondensatoren mit Differentialbriicke 


geschaltetem Verlustwiderstand, im anderen Zweig hingegen ein Kondensator 
mit parallelgeschaltetem Verlustwiderstand, kann die Briicke nur bei einer 
einzigen Frequenz abgeglichen werden (Bild 3). 
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Im folgenden wird vorausgesetzt, daB die Briicke bei der Grundfre- 
quenz abgeglichen wurde. Hine der zu vergleichenden Impedanzen hat Reihen-, 
die andere Parallelverluste. Der Verlustfaktor der beiden Impedanzen bei der 
Grundfrequenz ist 


D=tgo. 


Setzt man voraus, daB die Ohmschen Komponenten frequenzunab- 
hangig sind, die reaktiven Komponenten jedoch mit der Frequenz pro- 
portional (Induktivitaten) oder umgekehrt proportional (Kondensatoren) 
wachsen, ist die Ausgangsspannung, der n-ten Harmonischen mit der 
_ Amplitude k,U, 

i: 

pay saga 

pan EU, n 
n D) 


1 1 
DS i (n — — 
pi n 


Der Absolutwert der Ausgangsspannung ist 


n— — 


|U,| = =p 
(abel 


Bild 4 zeigt den fiir den Absolutwert der Ausgangsspannung charakteri- 
stischen Quotienten | U,|/k, Up. 


Wird ein phasenempfindlicher Indikator beniitzt, ist nur diejenige 
Ausgangsspannungskomponente von Bedeutung, die mit der Speisespannung 
in Phase liegt, in obigem Falle also 


1 1 
D + —||n — — 
= k,U, | D nm 
HC 2 
2 12 1\? 
D + —|+ jn —— 
| os D | n 
Diese Spannungskomponente erreicht bei 
1 
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ein Maximum. In diesem Falle wird 


= k,, Uy 


Bild 5 zeigt die auf die Harmonischen bezogene Ausgangsspannungs- 
komponente, die mit der Speisespannung in Phase liegt, in Abhangigkeit 
von dem auf die Grundharmonische bezogenen Verlustwinkel. 


Up) 
las 


Bild 4. Amplitude der Harmonischen in der Ausgangsspannung einer Differentialbriicke, 
die auf die Grundharmonische abgeglichen ist. 


001 Q02 O05 af a2 05 7 2 tg 


Bild 5. Komponente der Harmonischen, die mit der Speisespannung in Phase liegen, im Aus- 
gangssignal einer auf die Grundharmonische abgeglichenen Differentialbriicke 


Beim Messen von Induktivitaten kénnen die Harmonischen zufolge 
der Eigenkapazitat der Spulen sehr stéren. Die Briicke ist bei Grundfrequenz 
abgeglichen, wenn die effektive Induktivitat in beiden Briickenzweigen den 


EINE NEUE PHASENEMPFINDLICHE INDIK ATORSCHALTUNG 339 


gleichen Wert hat. Die wirksame Induktivitat einer Spule mit der Induktivitat 
L und der Streukapazitat C betragt 


ie L L 
effe —_ 3 — 9 
1 — w? LC rea (y 
worin Wo 
1 
oz = —— 
oglC 


die Eigenfrequenz der Spule ist. 
Es soll nun jener Extremfall untersucht werden, bei dem die eine Induk- 
tivitat vollkommen frei von Kapazitaten ist, wahrend im anderen untersuchten 


10 - 
Qo01 aor af f 


Bild 6. Amplitude der Harmonischen in der Ausgangsspannung einer Differentialbriicke, 
die auf die Grundharmonische abgeglichen ist 


Fall eine Streukapazitat C besteht, die man sich als zu der zu messenden 
Induktivitat parallelgeschaltet vorstellt. Die Kreisfrequenz der Grundwellen- 
komponente der Speisespannung sei @,. Wird die Bezeichnung 


Q=5 
ee 
eingefuhrt, ist die Abgleichbedingung 
L 
Ly = ——*—. 
pmele 1? 


Nach Abgleichen auf die Grundharmonische liefert die n-te Harmonische 


eine Ausgangsspannung 
_ Ip 270 — n?) 


he 
2 2+22(1—n?) 


Bild 6 zeigt die Briickenausgangsspannung als Funktion der Ordnungszahl 
der Harmonischen bei verschiedenen Eigenfrequenzen der Spule. 
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4. Beschreibung einer neuen Speise- und Indikatoreinheit [4] 


Im folgenden soll die neue schaltungstechnische Lésung einen MeB- 

briicke mit phasenempfindlichem Indikator beschrieben werden. Die Schalae 

ist bedeutend einfacher als die der tiblichen Lésungen, und auBerdem benétigt 
sie nur zwei Réhren mit je zwei Steuergittern (Bild 7). 

Die Réhren werden doppelt ausgeniitzt. Bei Pentagridréhren arbeitet 
das System drittes Gitter-Anode als Gegentaktoszillator, dessen Kreisindukti- 
vitat im Schwingkreis im Anodenkreis aus der Primarwicklung des Trans- 
formators besteht. Die symmetrischen Sekundarwicklungen dieses Transforma- 
tors bilden die beiden Zweige der Differentialbriicke, wahrend sich in den 


; +400V 


Bild 7. Schaltung der MeBbriicke mit phasenempfindlichem Indikator 


beiden anderen Zweigen die zu messende und die bekannte Impedanz befinden. 
Weicht der Wert der zu messenden Impedanz von demjenigen der Normal- 
impedanz ab, entsteht eine Fehlerspannung. Die GréSe der mit der Oszillator- 
spannung in Phase liegenden Fehlerspannungskomponente hangt von der 
GréBe des Unterschiedes der beiden Impedanzen, die Richtung der Fehler- 
spannung von der Richtung des Unterschiedes ab. Die Fehlerspannung wird 
an die beiden miteinander verbundenen ersten Steuergitter der Réhren gelegt. 
Da die Spannungen am ersten und dritten Gitter der beiden Réhren gleiche 
Frequenzen haben, wird die eine Komponente der multiplikativen Mischung 
Gleichstrom sein. Die Oszillatorspannungen an den dritten Gittern der beiden 
Réhren sind zueinander um 180° in Phase verschoben, die Gleichstromkompo- 
nente wird sich somit zum Ruhestrom der einen Réhre addieren, wahrend sie 
sich von dem der anderen Rohre abzieht. Schaltet man also zwischen die beiden 
Anoden ein Gleichstrominstrument, wird dieses nur den Anodenstromunter- 
schied messen, der wieder vom Unterschied der an die beiden ersten Réhren- 
gitter gelegten Wechselspannung, also vom Unterschied der beiden Impedanzen 
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abhangt. Das Instrument kann mithin direkt auf prozentuale Abweichung 
geeicht werden. Die gute Stabilitat der Anordnung wird durch den symmetri- 
schen Aufbau und durch die negative Riickkopplung am Kathodenwider- 
stand gesichert. Die Schaltung kann auch mit zwei Trioden-Hexoden aufgebaut 
werden. Die beiden Trioden bilden dann den Oszillator, der die Briicke speist. 
Die Schaltung ist der in Bild 3 gezeigten ahnlich. Die Verwendung von Trioden- 
Hexoden hat den Vorzug, daB der Generator und der Detektorteil in gewissem. 
Sinne voneinander unabhangig sind. 


5. Phasenempfindliche Anzeige [5], [6], [7] 


Man lege an das erste Gitter einer Réhre mit zwei Steuergittern die 
Spannung 


u, (t) = U,,coswt+ U,,sinat. 
Da am dritten Gitter der Réhre die Wechselspannung 
ug (t) = U, cos (wt + kx) 


liegt, ist die auf die Anode bezogene Steilheit des ersten Gitters zeitlich ver- 


Sy 


Bild 8. Zeitlicher Verlauf der Steilheit 


anderlich. Ohne Phasenverschiebung ist diese Steilheit eine eindeutige Funk- 
tion der Spannung des dritten Gitters (Bild 8). 
Der Anodenstrom ist 


ig (t) = S() 1 (1)- 
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Da S(t) eine ganze periodische Funktion ist, schreibt sich ihre FourreRsche 
Reihe 


S(t) = Sy + S, cos (wt + kz) + S_cos 2 (wt + km) + ary 


'S,cosn(wt+ kz). 


0 


ll Wa 


n 


Die Koeffizienten der Fourrerschen Reihe kénnen aus der Charakteri- 
stik S—U,, graphisch oder analytisch bestimmt werden. 

Setzt man die FourrERsche Reihe der Steilheit in die Formel des Anoden- 
stromes ein, erhadlt man 


i, (t) = (Uy,.cos ot + U,,sinot) S' S,cosn(ot+kz), 
n=0 


und nach entsprechender Umformung dieses Ausdruckes ]48t sich der Anoden- 
strom in Komponenten zerlegen, deren eine die Summen, der andere die 
Differenzen der verschiedenen Frequenzen enthalt. Das Kosinusglied der 
Differenzen bedeutet im Falle n = 1 Gleichstrom, dessen GréBe 


i,= 5 5:Uj- 008 ka. 


In der Schaltung sind die Phasen der an das dritte Gitter der beiden Réhren 
gelegten Spannungen entgegengesetzt 


k®) — 0 
k@) = 1, 


weshalb sich der Anodenstrom der einen Réhre verkleinern, der anderen Réhre 
vergréBern wird. 


Die beiden Gleichstromkomponenten sind also 


1 1 
I? — Sy i 1.008 0 = eee 


und 


1 1 
te Ct cos 7 = — Fy 


Ic? 


wahrend der Unterschied zwischen den beiden Anodenstrémen 


AI,=I9 — 19 = S,U,, 
betragt. 
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Da die Mischsteilheit 


ein im Katalog angegebener Réhrenkennwert ist, kann unter Benutzung 
desselben der Anodenstromunterschied 


LE ROIS is. 


leicht errechnet werden. 


6. Die Wirkung der Harmonischen 


Die Harmonischen in der Wechselspannung die an die verbundenen 
ersten Gitter der beiden Réhren gelegt wird, kénnen MeBfehler verursachen. 
Da die Briicke durch einen Transformator im Gegentakt-Anodenkreis gespeist 
wird, findet man in der Speisespannung der Briicke infolge der Gegentakt- 
schaltung beinahe ausschlieBlich Harmonische mit ungeraden Ordnungszahlen. 
Es soll nun untersucht werden, wie empfindlich sich der Indikator den Harmo- 
nischen gegeniiber verhalt. Hierbei wird vorausgesetzt, daB die Charakteristik 
der Mischréhren als bilinear angesehen werden kann, da8 also 


Tq = @ (U gy — Ugo) (u gg — Ugo). 
In diesem Falle schreibt sich die auf die Anode bezogene Steilheit des ersten 
Gitters 
Ss ELeoiey, (u 
O Uo1 


— Usp) . 


Die Vorspannung des dritten Gitters wird nun so hoch gemacht, daB der 
Arbeitspunkt in den unteren Knick (U4) der Charakteristik zu liegen kommt. 
Die Amplitude der Wechselspannung am dritten Gitter ist gleich der Vorspan- 
nung dieses Gitters. Diese Einstellung entspricht der Normaleinstellung von 
Mischréhren. In diesem Falle ist die Steilheit als Funktion der Zeit 


J 1 
a Ug) cos@t = Sm,xcos@t wenn — ey = Obs + i 
Z 
AU Beez Dt ee 
| | EL 
| wenn / 
O | fe) SOU SS GE 
2 = == 
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Die Fourrer-Koeffizienten der Steilheit sind 


+5 1 


Simca wenn n= 0 


bo 


5 
s, = "ms | cosmtcosnatd(wt) =, 
1 


Sis wenn n= 1 


— + ] 
(1) , wenn n= 2,4,..2k 
a 


ee 


0 wenn n = 3,5...2k+1 


Ist die Steuerspannung nicht verzerrt und die Charakteristik linear, ist 
das Indikatorsystem nur den geraden Harmonischen gegeniiber empfindlich. 
Da jedoch in Wirklichkeit das Oszillatorsignal verzerrt, und auch die Réhren- 
charakteristik nicht vollkommen linear ist, wird sich die Zeit-Funktion der 
Steilheit ebenfalls deformieren. Ein Spezialfall dieser Deformation besteht 
darin, daB an Stelle der halben Sinuswellen Dreieckimpulse entstehen. Der 
zeitliche Verlauf der Steilheit ist dann durch den Zusammenhang 


1 chee nm 
S() = Swax[ a+ 3 aoe 1 — cos 3 |cosmot| 
gegeben. In diesem Falle fehlen nur die 4., 8., .... Harmonischen. 


In einem anderen Grenzfall der Signaldeformation entstehen an Stelle 
der halben Sinuswellen Rechteckimpulse ; es gilt dann 


S (t) = Mare 


n=12n 


1 co 
Ss ~ 3a cone: 


wobei nur ungerade Harmonische auftreten. 

Weicht das Signal und die Charakteristik von der idealen Form ab, 
enthalt die Fourrersche Reihe der Steilheit auch ungerade Glieder, welche 
mit den ungeraden Harmonischen des Gegentaktoszillators eine Gleichstrom- 
komponente ergeben, die von der Phasenlage abhangt. 

Ist die relative Abweichung der zu messenden Impedanzen A, so tritt 
an den Ausgangsklemmen der mit U, gespeisten Briicke eine parallele Span- 
nungskomponente der Gréfe 


U 
U,.= A 
4 
auf, was einen Unterschied der Anodengleichstréme von 


AI, =S,Uy= 8,-2.4 


bedeutet. 
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-Enthalt die Speisespannung der Briicke eine n-te Harmonische mit der 
Amplitude k, U,, so ist — in dem in Punkt 2 behandelten unginstigen Fall — 
die parallele Ausgangsspannungskomponente der n-ten Harmonischen 

k,, Uy 


One ES : 
4 


Diese Spannung ruft einen Anodenstromunterschied von 


A Tan = SU e SE 


~~ hervor. Im ungiinstigsten Falle also (bei Addieren der Absolutwerte) erzeugen 


die Harmonischen einen Stromunterschied von 
eZ Uy, = 
Alans > mare > Sik 
n=2 n=2 
was einer Impedanzabweichung von 


co 
5 
H > 


n=2 


1 
entspricht. Dies ist somit der eréBte MeBfehler, den die Harmonischen erzeu- 
gen. 

7. MeBergebnisse 


Bild 9 zeigt die Schaltung eines fertigen Musters der oben beschriebenen 
VergleichsmeBbriicke, Bild 10 die AuBenansicht des Musterstiickes. 


Im weiteren sollen nun die mit der MeBbriicke erzielten MeBergebnisse 


erértert werden. 


ECG2 ; ECH 8 couse _+100¥ 


Bild 9. Schaltung der VergleichsmeBbriicke 
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Bild 11 veranschaulicht den Ausschlag des Instruments als Funktion 
der relativen Abweichung zwischen den zu vergleichenden Widerstanden. 

Der Spannungsabfall, der am Innenwiderstand des Differentialtrans- 
formators und des Oszillators entsteht, verringert beim Messen kleiner Wider- 


Bild 10. a) Vorderansicht der Vergleichsmefbriicke 


b) Innenansicht der Vergleichsmefbriicke 


stande die an die Briicke gelangende Spannung. Die Empfindlichkeit der 
Briicke ist deshalb bei derartigen Messungen etwas geringer. 
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Das Diagramm 12 zeigt den Zusammenhang zwischen der Empfindlich- 
keit und dem zu messenden Widerstand in einer Mefeinrichtung ohne Ver- 
starkerstufe. 

Es 1a8t sich feststellen, daB die Belastung durch den zu messenden 
Widerstand bei einer einfach ausgefihrten Briicke bei Widersténden von 


Ce ZT 
Ae 


Bild 11. Zusammenhang zwischen Instrumentausschlag und relativer Widerstandsabweichung 


10 kQ die Empfindlichkeit bereits bedeutend herabsetzt. Die Abhangigkeit 
der Empfindlichkeit von der Belastung kann durch Verringern der Uberset- 
zung des Differentialtransformators aufgehoben werden, sollen jedoch Wider- 


100 200 500 
10 20 50 0. 712) 


Bild 12. Zusammenhang zwischen Empfindlichkeit und Belastung 


stande in der GréS8enordnung von einigen Ohm gemessen werden, wobei die 
Briickenspannung stark abfallt, mu8 zum Ausgleich die Fehlerspannung durch 
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Z.wischenschalten einer Verstarkerstufe auf ein dem phasenempfindlichen 
Indikator entsprechendes Pegel gehoben werden. Mit Hilfe einer einzigen 
Verstarkertriode mit der Verstarkung A sinkt die zu messende kleinste 
Impedanz auf den A®-ten Teil des urspriinglichen Wertes, vorausgesetzt, 
daB die Ubersetzung des Differentialtransformators gleichzeitig auf den A-ten 


2 r/ 
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Bild 13. Zusammenhang zwischen Empfindlichkeit und Belastung bei Verwendung 
einer Verstarkerstufe 


Teil des urspriinglichen Wertes verringert wird, Bild 13 zeigt den Zusammen- 
hang zwischen Empfindlichkeit und Belastung bei 40facher Verstarkung. 
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Bid 14, Nullpunktverschiebung infolge 
Netzspannungsschwankungen 


Bild 15. Empfindlichkeit in Abhangigkeit 
von Netzspannungsschwankungen 


Die Nullpunktstabilitaét kann der an den Kathodenwiderstanden auf- 
tretenden negativen Riickkopplung und des s 


’ ymmetrischen Aufbaues wegen 
als sehr gut bezeichnet werden. 
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Bild 14 stellt die Nullpunktverschiebung dar, die infolge der Netz- 
spannungsschwankungen auftritt, auf die die gleichen Zeigerausschlige bewir- 
kenden relativen Widerstandsabweichungen bezogen. 

Die Netzspannungsschwankungen beeintrichtigen auch die Empfind- 
lichkeit, weil sich die Oszillatorspannung andert. Bild 15 veranschaulicht die 


ys torspannung 


_ Bild 16. MeBanordnung zur Untersuchung der Wirkung von Stérspannungen 


relative Impedanzabweichung, die den Endausschlag des Instrumentes bewirkt, 
in Abhangigkeit von der Netzspannungsschwankung. 

Es wurden weiterhin Messungen mit der Schaltung gem4B Bild 16 
durchgefiihrt, um den Einflu® der an der zu messenden Impedanz liegenden 
Stérspannungen auf die Arbeitsweise der Briicke festzustellen. 

Solange am ersten Steuergitter keine Ubersteuerung infolge der Stér- 
spannung auftritt, ist das Gerat praktisch unempfindlich gegen Stérspannun- 
gen, deren Frequenz von derjenigen der Mefifrequenz abweicht. 
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Bild 17. Wirkung der Stérspannung 


Im Diagramm des Bildes 17 sind diesbeziigliche MeBergebnisse auf- 


getragen. - 
Die Harmonischen des Oszillators und die nichtlineare Charakteristik 


der Mischréhre bewirken eine Reaktion des Indikators auch auf die Harmo- 
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nischen. Bild 18 veranschaulicht die Empfindlichkeit des Indikators fiir die 

Harmonischen bei 25 mV. Es fallt auf, daB® die vierte Harmonische praktisch 

wirkungslos ist. Untersucht man das Signal des Oszillators, so kann man fest- 

stellen, da® es sich der ungeraden Harmonischen wegen an beiden Scheiteln 
1 


5} 
St 


Of 


001 
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Bild 20. MeBfehler, hervorgerufen durch die Harmonischen 


zuspitzt und dem im Punkt 6 beschriebenen Dreiecksignal ahnlich gestaltet - 


Infolge dieser Signalform entsteht in der Steilheit keine vierte harmonische 
Komponente. 
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Bild 19 zeigt die Harmonischen im Ausgangssignal des Oszillators.. 
Die hier folgende Tabelle gibt das Verhaltnis der Mischsteilheit der einzelnen 
Harmonischen zur Mischsteilheit der Grundfrequenz, den Gehalt an Harmoni- 
schen und den fur die Wirkung der Harmonischen charakteristischen Wert 


Se 
S; 
an. 
n Sn/Sy kn (Sn/Sx) kn 
1 1,00 1,00 1,0 
2 0,08 0,01 0,8 - 10-3 
3 0,10 0,048 4,8 - 10 
4 0,04 0,0022 0,088 - 10-3 
5 0,12 0,013 1,56- 10-3 
6 0,02 0,00 0,00 
1 0,04 0,0026 0,104 - 10-3 
8 0,00 0,00 0,00 
9 0,04 0,00 0,00 


Im Bild 20 ist die far die Wirkung der Harmonischen charakteristische 
GréBe k,,S,/S, als Funktion der Ordnungszahlen der Harmonischen aufgetra- 
gen. 

In der beschriebenen Anordnung ist 


8 


eI k, = 7,25 + 10-8, 


1 


bv 


’ 
n=2 


der durch die Harmonischen verursachte prozentuale Fehler betragt also 
Aj, = 7,25 + 10-8 = 0,725 - 


Besteht ein Phasenunterschied zwischen der zu messenden und der 
Vergleichsimpedanz, so wird Nullindikation erzielt, wenn die Absolut werte 


Bild 21. Vektorbild der auf Absolutwerte abgeglichenen Briicke 
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genau abgeglichen werden. Aus den Vektoren in Bild 21 geht hervor, daB die 
Ausgangsspannung der Briicke in diesem Falle keine Komponente enthalt, 
die mit der Speisespannung der Briicke in Phase liegt, ist somit 


Al, =25,0,.=90. 
Gleichzeitig kann jedoch die senkrechte Spannungskomponente 


U @D 
U,,= Pry 


einen betrichtlichen Wert haben. Dieses Signal erzeugt Ubersteuerung, was 


0 20 40 60 80 $° 


Bild 22. Zusammenhang zwischen Phasenwinkel und Empfindlichkeit 


in der aus Bild 22 ersichtlichen Empfindlichkeitsverringerung in Erscheinung 
tritt. 


8. Anwendungsgebiete 


1. Grundlegend finden die phasenempfindlichen Indikatoreinheiten bei 
ImpedanzmeBbriicken Verwendung. Bei einem Teil der Briickenschaltungen 
geschieht der Abgleich komplexer Impedanzen derart, daB erst auf Absolut- 
wert, danach auf Phase abgeglichen wird. Ist die Vergleichsspannung des 
phasenempfindlichen Indikators bei dieser Art von Briicken wahrend der 
ersten Messung in Phase mit der Speisespannung der Briicke, bei der zweiten 
Messung dagegen um 90° verschoben, 1a8t sich der Abgleich an Stelle sukzessi- 
ver Annaherung direkt durchfihren. 

2. Briicken mit phasenempfindlichem Indikator kénnen bei der fabriks- 
maBigen Herstellung nachrichtentechnischer Bauelemente vorteilhaft ange- 
wendet werden, besonders zur Kontrolle der Toleranzen von Widerstinden 
und Kondensatoren. 

3. Die MeBbriicke mit phasenempfindlichem Indikator kann zur stan- 
digen Fabrikationskontrolle lackierter Kupferdrahte beniitzt werden [8]. Nach 
dem Patent Ingenieur H. Prunxu [9] lauft der lackierte Draht, dessen 
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Kupferleiter geerdet ist, durch ein QuecksilbergefaB das mit dem Draht 
einen Kondensator bildet, wobei die Lackschicht als Dielektrikum wirkt. Die 
in geeigneter Weise ausgebildete KapazitatsmeBbriicke wird auf die durch- 
schnittliche Lackdicke abgeglichen. Jede Abweichung der Lackschichtdicke 
vom Durchschnittswert lést dann einen Ausschlag des Instrumentes aus. 
Zwischen den Anderungen der Schichtdicke und der Kapazitat wurde bei den 
ersten Versuchen ein nichtlinearer Zusammenhang festgestellt, der jedoch 
guerst durch Verwendung einer Kristalldiode als Nebenschlu8 spater aber 
durch entsprechende Ausbildung der Polschuhe des Drehspulinstrumentes 
linearisiert werden konnte. 

4. Auf Grund einer Idee des Ingenieurs L. T6TH wurden Versuche mit 
einem batteriegespeisten kleindimensionierten Gerit zur Ermittlung des 
Stempeldruckes in Kohlenbergwerken angestellt. Zur Messung wurde ein 
MeBkondensator mit Gummidielektrikum benitzt. Durch diese Messung 
konnten bisher unbekannte Druckwerte festgestellt werden. ) 

5. Nach einer Idee des Ingenieurs T. ALKER [10] lie sich mit dem beschrie- 
benen Gerdt die Tiefe der geharteten Schicht von Stahlen feststellen : In 
die beiden Zweige der Differentialbriicke wurden induktive MeBkipfe einge- 
setzt. Der eine MeBkopf nahm das Werkstiick auf, dessen Hartungstiefe 
bekannt war, wahrend auf den anderen MeSkopf das Priifstiick aufgesetzt 
wurde. Die Induktivitat der MeBképfe hangt von der Harteschichttiefe ab, 
so daB auch diese Messung mit Hilfe des phasenempfindlichen Indikators 
durchgefiihrt werden konnte. 

MeBbriicken mit phasenempfindlichem Indikator kénnen bei der Lésung 
yieler technischer Probleme gut benutzt werden und sind besonders dort 
vorteilhaft, wo die zu lésenden Aufgaben auf die Messung von Kapazitaten 
oder gréBeren Widerstanden zuruckgefihrt werden kann. 
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Zusammenfassung 


Bei vielen wissenschaftlichen und technischen Problemen ist es ndtig, die GréBe und 
Richtung der verhiltnismaBig kleinen Anderung irgendeiner Impedanz festzustellen. Im 
allgemeinen benutzt man zu diesem Zweck WechselstrommeBbriicken mit phasenempfindlichem 
Indikator. Die Arbeit beschreibt eine neuartige Schaltung fiir die Speise- und Indikator- 
einheit und gibt einige Anwendungsbeispiele. 
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1. Production of enamel-covered wire 


In the machines and instruments of the electrical industry and of the 
electronics the insulation of the wires used may be different : one or more 
layers of textile, enamel cover, or a combination of the two. The advantages 
of enamelled wire against other insulated wires are its low production cost, 
relatively high heat resistance and its small space requirement. 

The basic material of the enamelled wire is generally copper, sometimes 
aluminium. In an apparatus, at a given power and temperature, the cross- 
section of the aluminium wire must be greater because of its lower conductiv- 
ity. If the insulation can stand higher temperature then aluminium wire of 
the same cross-section as copper can be used and this could be important in 
the question of price and weight reduction. 

Some time before, solely oleo-resinous enamel was used for enamelling. 
The fundamental material for this is the China wood oil. The development in 
the production of the synthetic materials has brought about important changes 
also in the enamels used for enamelling wires. Several different enamels of 
excellent quality have been developed, which could more or less be included 
in one of the following groups of plastics : 

polyurethane 

polyvinylformale 

polyesther 

copolymer of vinylchloride vinylacetate 

carbamide resin 

polyamide (perlon) 

These enamels have a much higher resistance against heat and chemicals, 
have higher mechanical strength etc., than those of oleo resinous ones. There 
is even now such an enamel which melts at the temperature of soldering and 
desoxidates the wire. This is especially advantageous in a multiterminal 
transformer, because the wire can be soldered to its place without preliminary 


cleaning. 
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The requirements for enamelled wires are different, depending on whether 
they are used in electrical or electronic engineering. The high number of 
turns common in electronic equipments demand that the thickness of the 
layer does not exceed the standard specifications. The continuity of the 
layer is also important, whereas by great dielectric strength is not, because 
of the small turn-to-turn voltage. Low dielectric loss is required of wires used 
in high frequency coils. 

In the machines and equipments of electrical engineering more important 
points are high dielectric and mechanical strength of the insulation and its 
resistance against heat and oil. The trend towards the specified load emphasize 
these factors all the more. 

The enamel layer must stick fast to the wire, this can be achieved by 
producing clean surfaces by using the proper enamel quality. The enamel layer 
must be flexible enough not to crack or split while spooling, because of the 
pulling and bending stress. The layer even under impregnation must be so 
hard that no short circuit should occur when two crossing turns cut into each 
other. 

The method of enamelling is based on putting more and more thin 
enamel layers — the thickness of which is of the order of microns — on the 
wire, burning it in on an electric heated oven, then winding the finished enamel- 
led wire on spools. The use of more than one layer is necessary, because the 
dissolvent must easily evaporate from the enamel-film. Should the layer be 
too thick, the dissolvent would evaporate with a burst and damage the smooth 
surface of the layer. 

The proper temperature for burning in is secured by automatic temper- 
ature regulating of the oven. The colour and glossiness of the enamel layer 
indicates the proper stage of the stoving. 

The thickness of the enamel layer is regulated by standard specifications. 
It is in the interest of the producing firm to stop the production of wires not 
fulfilling these specifications or to modify the thickness of the layer properly, 
even in the case of small deviation. 

Wires thicker than 0,5 mm move relatively slowly between the stove 
and the winding spool and there is time enough for someone to check the 
diameter with a micrometer by slowly walking along the moving wire. 

On the other hand, wires 0,05—0,5 mm in diameter are made by other 
types of machines, the speed of which are much greater. The greatest trouble 
in the micrometer measuring method is the necessity of again and again 
tearing the wire for measuring, until the enamelling is properly set. Moreover, 
the layer-thickness once adjusted does not remain constant either. As at a 
spool only once is there an opportunity to check the diameter — when taking 


the spool off — the thickness of the enamel layer cannot always be kept within 
the prescribed limits. 
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2. Continuous measuring of the layer-thickness 


Continuous measuring can be achieved by capacitive method PORTE 
the wire passes through a vessel filled with mercury, a cylindrical condenser 


Mercury Faeel c Copper 


Fig. 1. Measurement of layer-thickness by capacitive method 


is formed, one electrode of which is the wire itself, the dielectric is the insulat- 
ing layer, while the other electrode is the mercury (Fig. 1). 
The capacitance can be calculated from the equation 


Fee ss eG (1) 


(impo 
1 


where ¢e is the dielectric constant of the enamel, I the length of wire in cm 
merged in the mercury, d, the diameter of the base wire and d, the diameter 
of the insulated wire — these last two quantities should be substituted as the 
same unit otherwise any arbitrary measuring can be used. 

The dielectric constants of the enamels mentioned above are between 
2,8 and 4,0 at 800 c/s. The dielectric constant of an unknown enamel can be 
determined by way of the mercury filled measuring vessel, if we know the 
geometry of the wire. 

Since the Standards in most cases contain the diameter of the copper 
and the layer-thickness, it is advisable to eliminate d, from the equation and 
have the layer-thickness (v) as the new variable. So — since d, — dj = 2v — 


eee acts [pF] 

v 

In| 1 + — 

dy 
For the conversion of Standards [2]—[5] or rather to make quick estima- 
tions, the equation above is not easy to handle. If we consider the fact that — 
according to most of the Standards — between 0,05 and 0,5 mm nominal 
diameters the minimum diameter growth is not less than 5% of the copper 
diameter, the maximum diameter growth is not bigger than 50%, we get a good 


approximation with the relation : 
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i a 7 : 
PRUE RH TM et Na [pF]. (2) 
v 


At the mentioned limits the value of the amount in brackets is 20,48 
and 2,48; the absolute value of the deviation from the exact relation never 
exceeds 0,015 (see Appendix). So the capacitance as function of the layer-thick- 
ness almost hyperbolically varies. 

The Standard for enamelled wires gives the mean layer-thickness for 
every copper diameter and that tolerance range, inside which the thickness 
of the layer is acceptable. From these data one can calculate the minimum, 
mean and maximum capacitance belonging to the maximum, mean and mini- 
mum layer-thickness. The character of the relation between the copper diame- 
ter and the capacitances is shown in Fig. 2. One must note the heavy fluctua- 


Fig. 2. Relation between the nominal diameter of standard wire and its layer capacity 


tion of the curves and besides, the two capacitance differences : C,—C,, and 
‘Cy—C, are not equal. The reason of the fluctuation is due to the fact that 
the Standard supposes mechanical measuring device (micrometer) and, 
therefore, it prescribes the maximum, mean and minimum thickness of the 
layer — often arbitrarily — in whole numbers of microns. The discrepancy 
of the two capacitance difference, on the other hand, arises, because the thick- 
ness tolerance is symmetrical to the centre, but the capacitance is not a linear 
function of the thickness. In the capacitive layer-thickness measuring devices 
now in use, it was easy to realize the asymmetrical measuring range, but in- 
stead of the zig-zagging tolerance limit, a narrower limit had to be used (dash- 
and-dot lines in Fig. 2) since it would have been much more complicated to 
set in the exact limits on the apparatus. 

From equation (2) one can calculate the capacitance variation correspond- 
ing with the minimum and maximum layer-thickness : 


AC, = Cy — C, = 0,551] 2 — a (3) 


2Um 20% 
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pe ES, 
—_ et 1 
AC, =C, = Cn = 0,551] emacs |: (4) 
VE y Um 


where Cy is the maximum capacitance belonging to Um, the minimum layer- 
thickness ; C; is the mean capacitance belonging to v,, the mean layer-thick- 
ness; and C,, is the minimum capacitance belonging to vy, the maximum layer- 
thickness. The whole capacitance variation is 


46, + 46, = 0.55 61 [ 8 é (5) 
UV, 20m 
and the asymmetry 
i agieg 
AG, ema | (6) 
ena s eens 


If the thickness tolerance is symmetrical, and so 
Uy — UE = VE Ym (7) 


we get a simple relation for the asymmetry : 


git! eae (8) 


If the measuring is made according to the capacitive method, it would 
be preferable to have the value of (5) and (8) constant. The latter is fulfilled 
in the well-known Standards but the first condition can be realized only if, 


according to: 


dy al ee (9) 


a Ue 2vy 20m \ Um 


is constant. According to the Standards, 


Zi or as it comes from above, 
22 UM Uk 
though, along with the increase of diameter the relative layer-thickness decreas- 
es. That is the reason why it is necessary in the already working instruments 
to use a special setting for the adjustment of the mean capacitance — that 


is : of the nominal wire diameter. 


7 Periodica Polytechnica El II/4. 
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3. Constructing the ideal Standard 


It is advisable to choose the ratio of the layer-thickness and the diameter, 
so that the electric load of the wire of different diameters 


2U 
2v 
ayn [1+] 


a 


E= 


should be constant. With this solution it could easily be secured to have the 
outstanding layer-thicknesses eliminated, that is the over- and under-rated 
insulations. 

The Standards [2]—[5] presently in use are — as mentioned before — 
accomodated to possibilities given by the mechanical measuring devices 
(micrometer). The layer-thickness is given in whole numbers of microns. It is 
to be noted that the relative layer-thickness decreases with the increase of 
the nominal diameter. As the steps of the layer-thickness are chosen from 
different aspects (easy to memorize, symmetrical tolerance range, etc.), the 
relative layer-thickness as function of the diameter, shows strong fluctuations 
instead of the desired steady character. 

In order to realize a simple plant measuring method, it is more advisable 
to have such a Standard, in whichthe steps of diameter values are more or less 
smoothed to the capacitance curve. When forming the standard let us make 
the following conditions : 


1. The difference of the maximum and minimum capacitance should 
be constant. 

2. The ratio of the deviations from the mean capacitance should be 
constant. 

3. The relative layer-thickness should decrease with the diameter, so 
that the electrical load should remain nearly constant. 


All three conditions can be fulfilled only if we abandon the symmetrical 
thickness tolerance. The minimum, mean and maximum layer-thickness can 
be determined from the following equations : 

from the Ist condition : 


== ke 
Deny ee eG (10) 


from the 2nd condition : 


dq, da, =h| = 


2Un 2, 


(11) 


20, Jy 
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From the combination of (10) and (11) 


d, = d, aT ky (12) 
20yj 20, .1+hk, 


and 
d, d, ky kp 


= (13) 
20, 2%, +k 


d 
The value of 5 1_ can be given, either by taking the 3rd condition, or by the 
VK 


conventional Standards so far used. 

Fig. 3 shows the relation between capacitance and layer-thickness, 
according to the Hungarian Standard MSZ 1569. In the lower part of the 
figure a stepped diagram shows the layer-thickness as function of the diameter. 


CpF 
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Fig. 3. Relation between layer thickness and capacity according to Hungarian” Standard 
MSZ 1569 


The curves in Fig. 4 were plotted in the following way : we prepared 


the function “ osculatory to the middle curve in Fig. 3 ; then we determined 


the value of k, and k, from the positions of the exterior curves 


oe 5,49 
0,55 el 
ky ss AC, = 1,62 ' 
AC, 


7* 


362 A. AMBROZY, H. PRUNKL and K, TARNAY 


The values of vy and v,, from (12) and (13), rounded to whole microns, 
give the exterior curves of Fig. 4 and the stepped diagram, in which we can 


ss Qt G2 03 04 Q5 dmm 
2y% 
40 
20 
a1 Q2 a3 Q4 ~ 05 aymm 


Fig. 4. Relation between layer-thickness and capacity according to ideal Standard 


see, that we can design a standard very close to the ideal one, if we introduce 
a small change in the tolerance range. 


4, The principle structure of the measuring apparatus 


The main outlines of the measuring apparatus are shown in Fig. 5 
The measuring is based on bridge-method : an oscillator with the approximate 
frequency of 1 Ke/s feeds the primary windings of a differential transformer. 


Measured wire 


Fig. 5. Schematic block of measuring apparatus 


The mercury of the measuring vessel is connected to one end of the centre- 
tapped secondary coil, that is the outer plate of the capacitance (C,) to be 
measured. The other plate of the capacitance is the copper wire, which is 
grounded. A variable capacitance is also connected to these same points. 
The other end of the secondary coil is connected to a fixed capacitance. 
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The dial of the variable capacitance is calibrated to the values of the 
mean layer-thickness (v;), so that when measuring a wire of the diameter d,, 
we balance the bridge by setting the value of the capacity 


Co 60, a 0pbel te =f 0,48| [pF] (14) 
Ry 


which comes from equation (2). 
The capacity of the measuring vessel, if the diameter of the merging wire 
is d,, and its layer-thickness v, 


CG, 2055 ei ee 0,48| (15) 
v 


The output voltage of the differential-bridge between the centre tap 
of the transformer and the ground is, by using equation (14) : 


AC 
U, =U i C+, | 1 U C,—G, = Boks 2 oe (16) 
Tee (Pas Co teCn [ie oe CSC 2 ae 
C 


Finally, equation (3) gives the connection between the capacity variation 
AC, and the deviation of the layer-thickness. 

The phase of the voltage may be the same, or the reverse in respect to 
the bridge supply voltage, depending on whether the layer thickness is greater 


Fig. 6. The output voltage of the bridge plotted against the capacity variation 


or smaller than the mean value. The relation between the voltage and the 
capacity variation is plotted in Fig. 6. 

We can indicate the output voltage of the bridge by a vacuum tube volt- 
meter. The disadvantage of using a common type vacuum tube voltmeter is, 
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that it only indicates the amplitude of the voltage (Fig. 7.), but does not give 
information on its phase. That is the reason why we have to use a phase-sensi- 
tive voltmeter. 

We can get the reference voltage needed for the phase sensitive voltmeter 
from the same source, where we get the bridge supply voltage from. In this 


-10 -5 5 bCy 10 
4 


Fig. 7. The absolute value of the output voltage of the bridge, plotted against the capacity 
variation 


way the polarity of the output DC voltage corresponds to the direction of the 
layer-thickness deviation, while its value is proportional to the greatness of the 
deviation. The indicating miliammeter of the phase sensitive voltmeter is 
fitted with specially shaped pole-shoes (Fig. 8a) in order to equalize the non- 


Fig. 8. Miliammeter with specially shaped pole-shoes a) outline, b) sensitivity character 


linearity existing between the layer-thickness and the bridge voltage. So it 
has an asymmetrical sensitivity character (Fig. 8b). 

In this way it can be achieved that the indication of the meter is an almost 
linear function of the deviation of the layer-thickness from its mean value. 


d. Plant measuring equipments 


The block diagram of the first plant measuring equipment, of this kind, 
is shown in Fig, 9. 

A 6AU6 miniature pentode generates an A.C. voltage of 800 c/s, as a 
phase-shift oscillator. This is amplified by a 6AQ5 power pentode tube. A trans- 
former, in its plate circuit, supplies the differential transformer and gives the 


reference voltage for the phase detector using a 6AL5 twin diode (Fig. 10) [6]. 
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In one arm of the differential bridge is a fixed capacitance, in the other a 
variable one, the sum of the capacity to be measured and the other adjusting 
capacitances together have exactly the same value, as the former fixed capa- 
citance. In place of the capacity to be measured any of the 8 measuring chan- 


Oscillator 


Amplitier 


Fig. 9. Block-scheme of plant measuring apparatus 


nels of the 8 simultaneous enamelling channels can be switched, one after the 
other. The output voltage of the bridge is directly connected to the input 
terminals of a two-stage amplifier in the case of double layer. If, however, 
the wire is covered by a simple layer, a voltage divider is used. The amplified 
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Fig. 10. Phase-sensitive detector 


error signal goes on to the phase detector, which circuit can be seen in Fig. 10. 
In this version the linearity of the scale is not secured with specially shaped 
pole-shoes as mentioned above, but by a shunt crystal diode. The supply 
voltage needed for the instrument is taken from a mains-rectifier stage using 
a 6X4 miniature vacuum tube. In the Enamelled Wire Works (Albertfalva) 
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instruments of this kind have been already used for several years and they 
made it possible to produce more enamelled wire of better quality. 

By using the performance data and by introducing a new phase sensitive 
circuit stage [7]—[8], a new type of measuring instrument has been made for 


ECH81 Overall sersitivity ECH&1 


100 V DC. 
Supply 


Fig. 12. Photo of new cireuited plant measuring apparatus 
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export purposes, which contains only two vacuum tubes of the same type, 
besides the rectifier tube. The low number of tubes reduces the price of the 
instrument and greatly increases its reliability of service. The circuit diagram 
is shown in Fig. 11, while its photo can be seen in Fig. 12. Its working prin- 
ciple is as follows : 

The oscillating systems of two mixer tubes (ECH 81) work in a push- 
pull oscillator circuit, so the mixer grids receive voltages of opposite phases. 
The transformer in the plate circuit of the oscillator supplies the previously 
discussed differential bridge. The output voltage of this controls the joined 
first grids of the mixer tubes. Multiple mixing takes place and the frequency 
difference results in direct current, which will increase the DC current in one 
of the tubes and a decrease in the other, because of the opposite phases of 
the voltages connected to the oscillator grids. A capacity deviation in the 
opposite direction will result in a plate current difference in the opposite 
direction. The plate circuits of both tubes form a DC vacuum tube voltmeter, 
a midposition miliammeter with specially shaped pole-shoes is connected to: 
these points to indicate the unbalance. 


Appendix 


It is to be proved that equation (2) is a good approximation of equation 
(1) so 


: FY : + 0,48 
In (1 + x) % 
if 0,05 <x <.0,5. Starting from equation (1) we can get equation (2) in the 
following way : 
Let us expand in Taylor series the function In (1 + x) at the x = 0 place = 


x2 = x8 
In(l1t+«)=*x-——+—>7" 
n(1 + x) ae 
So 
i 3 1 i} : 
lw! 22 eae DB “224 5e—..1|| 
2, 3 2 3 4 
By a binomial series expansion it can be written : 
2, 
S ws pte(—ferge--|+ 
In (1 + x) x 2 3 4 
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Let us take only the first two members of the series in [ ], and to continue 
the neglection, let us regard constant as the so far unknown multiplier of x/2. 
In this way : 

A 

a" 


; 1 
1 ey 1 ! A x | re 
In (1 + x) a 2 x 
The error of the approximative formulae is given by the 
| a a | i] 
h(x) =—+ 
x 2 In (1 + x) 


error function (Fig. 13). 


hfx} 


Fig. 13. Error function 


The two functions smooth the best to each other, when the error integral 


H (x15) = [ h(x) dx 


xy 


taken for the whole range is zero. This condition gives the zero place of the 
error, and the A multiplier, respectively 


H (x2) =f (s dx i ue ——() 
Ce 8 Mie In (1 + x) 


To solve the first two integrals is not difficult, but the third one can only be 
calculated numerically. This latter one is the logarithmic-integral function, 


which can be found in tables of higher functions [9] with the symbol : li (1+-x). 
So 


HT (erg Xo) =n a) | s (x, — x,) — [i(1 4+ x,) —li(l1+4%,)]=0 


x4 
and at the given boundaries : 
x, = 0,05 li (1,05) = —2,3935 


%, = 0,5 (1,5) = 0,1251 
In 10° 2" 23026 
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So it comes that: 


0,1251 + 2,3935 — 2,3026 
0,45 


Aree == 0,96. 


The error is the greatest at the boundaries : 

h(x,) = — 0,01, 

h(x) = + 0,015. 
There is no error if x is equal to the root of the following transcendental 
equation 


boa 1 
+= =0 
x 2 In (1 + x) 
Here x = 0,3. 


Summary 


Both the heavy-current and weak-current industry put up strong demands against 
the enamelled wires. The thickness of the enamel layer cannot be continuously checked, in 
production, with mechanical measuring devices. The paper discusses a new method and a 
plant measuring apparatus, which measures the layer-thickness non-destructively in a capa- 
citive way. In connection with the new method a few questions of standardization are also 
discussed in this article. 
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INDUSTRIAL REVIEW--AUS DER INDUSTRIE 


THE ROLE PLAYED BY THE GANZ FACTORY IN THE EARLY 
ELECTRIFICATION OF EUROPE AND THE PRESENT ELECTRIFICATION 
OF THE BALKAN AND LEVANTINE COUNTRIES 


The Hidroelectric Power Plant of Ikizdere 


It was during the last decades of the 19th 
century that a world-wide interest in the 
new source of energy, namely electricity, 
was first aroused, and this interest has been 
steadily growing ever since. A bewildering 
variety of new machines, equipments, struc- 
tures and instruments has been brought forth 
with a view to utilizing, ‘‘taming” and trans- 
forming the new force, 1. e. electric current. 
In the course of this process, the dynamo, 
the telephone, the incandescent lamp were 
invented ; Déri, Blathy and Zipernowsky — 
all three of them engineers of the Ganz 
Factory — constructed the first transformer ; 
briefly the new energy — promising a pace 
of development never dreamed of before — 
began to conquer the world, and is still 
steadily gaining ground in the field of tele- 


communication and machine-tools. Electr- 
ical industry is enjoying a miraculous 
development. 


Endeavours to produce electrical energy 
in the cheapest possible way date from its 
very birth. It did not take long for attention 
to be directed towards hydraulic power which 
was only natural, for this form of energy — 
supplied gratuitously by nature — had been 
exploited by mankind since times immemo- 
rial. Side by side with the earlier directly- 
working turbines of fluor mills, paper mills, 
textile works and other industrial water 
turbines there appear the first hydroelectric 
power plants. 

Never since having brought out its first 
water turbine in 1866 has the Ganz Factory 
ceased to keep abreast with the progress of 
turbine construction. The initial period is 


characterized by the manufacture of Jonval 
turbines, Girard’s partial-admission turbines 
and, as a combination of these two types, the 
double-crown turbines. Such turbines, as 
produced by the Ganz Works, had a head up 
to 100 m and a capacity up to 700 h. p. 
The appearance of the Francis turbines meant 
a decisive turn in the process of development. 
Owing to comparatively high r. p. m. and 
Fink’s method of vane adjustment, this type 
proved particularly suitable for the driving 
of electric generators. It was in 1897 that 
the Ganz Works constructed its first turbine 
of this type. Francis turbines took a rapid 
development hereafter and, at the same time, 
also turbines of the Pelton type began to be 
increasingly used for the utilization of steep- 
er heads. 

As early as 1898 did the Ganz Factory 
supply the units required for the construc- 
tion of the hydroelectric plant at Jajce 
(Bosnia), i. e. eight 1000 h. p. and two 632 
h, p. Francis turbines with horizontal axles. 
Turbines made by Ganz could be found al- 
most at every point where electrification had 
gained ground, and it is not an exaggeration 
to say that the first really significant Europe- 
an hydroelectric power plants were the pro- 
ducts of the Ganz Works. It was likewise in 
1898 that Ganz delivered four 1400 h. p. 
turbines to the power station of the Etsch- 
werke at Bozen-Meran, four 1600 to 1800 
h, p. turbines to the power station of Tivoli 
which supplied the lighting of Rome, while 
— in 1899 — two 1200 h. p. turbines were 
manufactured for and delivered to the city 
of Innsbruck. After this, the old types — 
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Girard and the double-crown turbines — 
began to be obsolete, so that, for low and 
medium heads, nothing (or almost nothing) 
but Francis turbines had come to be em- 
ployed: these were single or twin turbines 
built in open shafts with a spiral waterway, 
or so-called boiler (twin) turbines. The North 
Italian plant at Morbegno, erected in 1900, 
was provided with three such 2000 h. p. 
Francis turbines. 

The four twin spiral Francis turbines, 
each with an output of 6000 h. p., manu- 
factured for the Dalmatian town of Manoilo- 
vac and the two 4000 h. p. turbines made 
for the Subiaco plant of the city of Rome 
in 1904, were really outstanding technical 
achievements at the outset of the present 
century. In addition to those we have men- 
tioned, the following of the several hundreds 
of power plants equipped with Francis 
turbines, deserve to be noted: that of Arci 
with two 5000 h. p. twin spiral turbines 
which supply the city of Rome with electric 
current, further the first construction of the 
Almissa power plant in Dalmatia (1911) 
which was provided with two twin spiral 
turbines of 20000 h. p. each. 

The Pelton turbines, developed along 
with the Francis type, are likewise signif- 
icant milestones on the path of progressing 
electrification in the history of the Ganz 
Works. Such turbines were made by Ganz 
for the Italian town of Spoleto in 1896 and 
then for the town of Foligno : they provided 
the current necessary for the lighting of these 
towns. A whole series of Pelton turbines were 
manufactured for the rolling mill of Kraina- 
Assling, among them a unit of 1738 h. p. 
for Jauerburg (head: 320 m). The four twin 
Pelton turbines of 1500 h. p. each (head: 
324 m) delivered to Port Madoc, England, in 
1904 merit special mention. The next deli- 
very to England was that of a 3000 h. p. 
unit in 1907 for the Conway Aluminium 
Works. 

Uncertain as the prospects of turbine 
production seemed to be in Hungary after 
World War I, the management of the Ganz 
Works decided not only to maintain but even 
to develop the manufacture of water-power 


machines. In taking this decision, the man- 
agement was encouraged, partly by the boom 
which it hoped would arise on account of 
the expected large-scale post-war construc- 
tions of hydroelectric power plants in the 
Balkan states, and partly by communica- 
tions and news concerning plans of a large- 
scale development of irrigation and drainage 
in Egypt. 

We must remember that, as far back as 
the last decade of the past century, it had 
become obvious that — reasonably — the 
manufacture of rotary pumps (performing 
the inverse of the work done by turbines) 
should be taken in hand, side by side with 
the continued production of turbines. It 
was, thus, at this time that the production 
of standard pumps, for mines and water- 
works in the first place, was begun. 

With a view to marketing both kinds of 
hydraulic machinery the management of the 
Ganz Works established commercial rela- 
tions with numerous countries around 1920: 
with Roumania, Bulgaria, 
Greece, Turkey, Italy, Spain, Poland, Egypt, 
India and various South-American countries. 


Yugoslavia, 


Outstanding among supplies made at this 
time were the construction of a hydroelectric 
plant at Patras (Greece) in 1924 with 1000 
h. p. Francis turbines (head: 115 m), as 
also two 30000 h. p. Francis turbines 
(head: 109 m) for the completion of the 
above-mentioned power plant of Almissa. No 
small credit may be derived from the fact 
that the turbines of the Almissa plant are 
still in full service and have never ceased to 
perform their work, without any general 
overhaul, to the complete satisfaction of the 
Dalmatian plant management. 

As proved by orders received and deliv- 
eries effected, between the two world wars, 
and after World War II, the Ganz Factory 
has always kept pace with the development 
of water turbines. 

It was especially the general industrial 
prosperity ensuing after the second world 
war which set significant new tasks to those 
engaged inthe production of water-power 
machines. Expanding industrialization re- 
quires more and more energy supply and 


prs 


ae 


gives rise to the erection of a steadily in- 
creasing number of hydro electric plants. The 
management of the Ganz works was well 
prepared for this onrush and is now in a 
position to manufacture most up-to-date 
equipments which are capable of satisfying 
all reasonable requirements. 

All water-power plants, as also many of 
the larger pumping stations made by Ganz 
during the last decades, were completely 
automatized. Some of them do not require 
any operating crews and can be controlled 
from a distant power plant. The latest pro- 
ducts are automatized to such an extent as 
to need only the pressing of a button for all 
necessary operations, including start and 
stoppage, to be performed, and it is only the 
load changes which have to be regulated from 
the control room. Noteworthy results have 
been achieved by the Ganz Works also in 
the development of speed-regulators. 

Since the termination of World War II, 
the Ganz Works has not only been able to 
meet Hungary's own demand of hydraulic 
machines but has also built up a considerable 
export market for its products in Roumania, 
Bulgaria, Poland, the Soviet Union, Turkey, 
India, Egypt, China and Yugoslavia. 

A new type of turbines, the 
Mignon” micro turbine, united with a gener- 
speed governor 


‘“*Ganz- 


ator and provided with a 
and voltage adjuster, was designed during 
the last years and is now manufactured in 
series. The whole structure constitutes a 
complete unit, its assembly and manipulation 
are quite simple, and is eminently suitable 
for the generation of current required by 
smaller communities, settlements and farms. 
Traditions almost a century old, together 
with theoretical and practical experiences, 
justify the claim of the Ganz Works to be 
ranked with the world’s most prominent 
manufacturers of water-power structures. 
It was a few years ago that the Iller 
Bankasi (an investment bank in Turkey) 
invited tenders for the construction of a 
hydroelectric power plant, that of Tkizdere. 
Among other tenders the most favourable 
and most complete was that submitted by 
Ganz which, on account of its undertaking 
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to construct and erect both the electric and 
hydraulic equipments, succeeded in obtain- 
ing the indent. 

“‘Ikizdere’’ means “‘twin brooks”: the 
name of this river derives from the fact that. 
it arises from the confluence of two equal 
branches. Strictly speaking, not just a single 
power plant but a whole system is being built 
at and in connection with Ikizdere, a place 
removed some 4.0 kilometres from the seashore. 
Situated in the northeastern corner of Turkey, 
this area in Asia Minor is flanked on two 
sides by the Soviet Union and the Black 
Sea, respectively. Supplying this region with 
electric current, the newly established system 
will satisfy an important requirement of cul- 
tural development in this part of the country. 
In addition to Ikizdere itself the said system 
embraces four littoral townships and their 
environments. The establishment in ques- 
tion must be regarded as a very signif- 
and its importance is the 
as no similar energy- 


icant one, 
more considerable, 
supplying system can be found within a radius 
of 500 kilometres. 

The technical task involved merits special 
mention. The water of the river, for instance, 
is conducted to the plant in two sections 
through a tunnel which, blasted in the gra- 
nite rock, has a length of 3% kilometres. 
Arrived at the plant, it is by means of a 
steel pressure tube of 300 m length and 2 m 
diameter that the water is dashed upon the 
turbine. 

The most important part of the machin- 
ery is a hydroelectric equipment consisting 
of three units. The water turbines are of the 
Francis type, with vertical axles and scroll 
cases which, utilizing a net head of 160 m, 
consume 4 m® per second of water and have, 
each, an output of 7500 h, p. with 750 r. p. m. 
Highest efficiency : 89 per cent. A like- 
wise vertically-axled three-phase synchroniz- 
ing generator is directly coupled with the 
turbine. Technical data of the generator: 
6500 kVA permanent output ; 6300 V volt- 
age; 50~; 96% efficiency. 

Regulation in accordance with varying 
operative conditions is effected by means of 
an automatic, oil-pressure speed-regulator, 
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system Ganz. A constituent part of the 
hydraulic equipment is a spherical valve 
before the turbine: it has a diameter of 1000 
mm and is likewise automatically governed 
by oil pressure. 


The electrical outfit includes the com- 
plete equipment and the instruments of the 
engine room, the transformer station and 
the switchgear in the open air, the trans- 
aission line and the substations. The current 
generated by the power plant is to be con- 
‘ducted by a 40 km long line to Iyidere on 
the seashore, therefrom to the town of Rize 
eastwards along the coast, as also to Arakli 
and Trabzon (and its surroundings), towns 


to the west. A step-down transformer station 
will be erected in each of these towns. 

The entire hydraulic and electric outfit 
of the power plant, including substations 
and transmission lines, is manufactured in 
and is to be shipped from Hungary. 

The construction at issue is a striking 
proof of the fact that, true to its traditions 
the Ganz Works has succeeded in maint- 
aining the good quality of its products, is 
capable of speedy performance, and is not 
inferior to even the most renowned and 
most prominent manufacturing establish- 
ments of the world. 

Compiled by 
Erwin Kopta 


THE EFFECT OF NOISE SOURCES ON THE PERFORMANCE 
OF A MULTI-CHANNEL FM RADIO LINK 


G. BoGnAr 


Introduction 


Microwave radio relay lines are designed 
with a most particular view to the fidelity 
of the signal transmission between the input 
and the output of the system. One of the 
most important characteristics of a multi- 
channel telephone system operated on a radio 
link is the total noise power appearing within 
one channel in the output signal. The length 
of the radio relay line is limited by the 
highest noise power which can be allowed in 
a telephone channel. Disturbances affecting 
multi-channel transmission are of ‘‘noise” 
character and can be divided, with a view 
to their origin, into three fundamentally 
different groups 

a) thermal noise, 

b) intermodulation between speech chan- 
nels, 

c) interference due to external sources. 

It can be shown that FM-systems are 


affected by two prevailing types of inter- 
modulation. These are the non-linear ampli- 
tude/frequency-deviation characteristic and 
the fluctuations of the group delay time. 
A characteristic feature of FM-systems con- 
sists in the fact that no distortion is caused 
by non-linear characteristics of tubes or non- 
uniform frequency response of the ampli- 
fiers. On the other hand, non-linear phase 
characteristics of the amplifiers, reflections 
in the antenna feed line and multi-path 
signal propagation cause time-delay distor- 
tions which may become critical if a number 
of relay stations is inserted. 

Multi-channel telephone systems, to eli- 
minate disturbances due to intermodulation, 
have to comply with high linearity require- 
ments. However, phase variations within the 
relatively narrow band of the individual 
speech channels can be tolerated. 


Sources of intermodulation noise 


The mostimportant non-linear distortions 
encountered in FDM—F™M multi-channel sys- 
tems can be classified, with respect to their 
origin, as follows : 


1. Distortion of the basic-band equip- 
ment due to non-linear characteristics of the 
amplifiers. 

2. Distortion of the modem equipment 
due to its non-linear amplitude and phase 
characteristics. 

3. Distortions of the high-frequency trans- 
mission path, which are 


a) distortions due to the finite band 
width, 

b) non-linear IF and RF phase charac- 
teristics, 

c) reflections in the antenna feed line, 

d) selective fading. 

With FM signal transmission two large 
groups of non-linear distortions can be dist- 
inguished, i.e. the static and dynamic distor- 
tions. Coefficients of static distortions are 
independent of the modulating frequency, 
while those of dynamic distortions increase 
with the same. 


Addition of intermodulation noises and their distribution aiong the hypothetical reference circuit 


When considering addition problems of 
intermodulation noise powers it is advisable 
to separate the examination of the problem 
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of addition between modem sections from 
that within a single modem section, because 
between the modem sections the phase 


3/0 


characteristic in the basie band is a prevail- 
ing factor. Noises due to crosstalk between 
modem sections can be assumed to have a 
power-summation law, the phases of the 
basis-band singals being related at random. 

Within a single modem section crosstalk 
disturhances due to amplitude and phase 
distortions are power-additive. Within one 
station the different sources of intermodul- 
ation noise can be considered as independent 
and for this reason, as a first approximation, 
an addition of the distortions due to different 
sources can be assumed. 


If considering the result of non-linear 
products of identical type produced within 
the same modem section no sitaple law of 
addition can be found. 

When designing for optimum frequency 
deviation the thermal noise and the inter- 
modulation noise are equal. Actually, the 
intermodulation noise increases with the 
frequency deviation, while the thermal noise 
decreases. Accordingly, the total psopho- 
metric noise in the output signal of the 
hypothetical reference circuit can be distri- 
buted as follows : 


10 000 picowatts 


-— Total noise == 


250 picowatts 
Noise of the basic- 
band equipment 


The ground noise is composed of thermal 
noise and radio frequency interferences in- 
cluding a fading reserve of 5 decibels. Inter- 
modulation noise power is to be determined 
by psophometric measurement over a band 
width of 3,1 ke/s. With the allowable inter- 
modulation noise power distributed among 
the nine modem sections, the noise power 
which may be generated within a single 
modem section is 420 picowatts. The follow- 


ing calculations will be based on unweighted 


wa — 210 picowatts 


105 picowatts 
amplitude nonlinearity 


52,5 picowatts 52,5 picowatts 


modulator demodulator 


According to the above distribution plan 
one half of the noise power should be due to 
amplitude distortion and the other half to 
phase distortion. Modulator and demodulator 
are supposed to have equal shares in these 
distortions. 


With the 630 picowatts, calculates for 


7500 picowatts 
High-frequency 


circuit noise 


3750 picowatts 
Ground noise 


3750 picowatts 
Intermodulation 
noise 


noise over a band-width of 4 ke/s and, 
to compensate for this, the admitted level 
of the intermodulation noise power will be 
increased by about 3 decibels to 840 pW. 
One quarter of this noise power will be as- 
signed to the modem equipment and three 
quarters to five transmitters and receivers. 
In this way the admissible noise power 
attributed to one*-modem equipment is 210 
picowatts which is divided as tollows : 


105 picowatts 
phase nonlinearity 


52,5 picowatts 
modulator 


52,2 picowatts 
demodulator 


the five transmitter-receiver pairs, the noise 
power which can be due to each pair is 126 
picowatts. As a first approximation, this 
noise power should be divided equally among 
the noises caused by feed line reflexions, 
non-linear IF-phase-characteristics and non- 
linear RF-phase-characteristics ; 


ale 


Sit 


126 picowatts 


a = 
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picowatts A2 


42 


Reflexions in feed line 


Distortions due to non-linear phase cha- 
racteristics in RF-sections include the phase 
distortions of the RF-filters and the disturb- 
ing phase modulation due to amplitude 
modulation in the travelling-wave tube. 


Relation between electrical characteristics 


and limit values of notse power 


When examining the effect of different 
noise sources on the transmission charac- 
teristics of the system, a method of measure- 
ment should be assigned to each source of 
noise. Amplitude nonlinearity are character- 
ized by harmonic distortion measurements, 
phase nonlinearity by measuring the fluc- 
tuation of the [IF and RF time-delay, and 
the reflexions in the transmission line by 
standing wave ratio measurements. In ad- 
dition, a measurement indicating the amount 
AM-signal 


amplification in the travelling wave tube is 


of phase distortion saused by 


also required. 

Noises due to the non-linear amplitude 
characteristics are interrelated with measured 
value of linearity through the distortion 
coefficient. With quadratic characteristies 
supposed as first approximation the following 


relation can be established between the noise 


power and the inverse of the distortion 
factor: 

=] 4 P? 

ie Y d, Ng 


where 


F is the amplitude of the test signal, 
H, is the amplitude of the second har- 


monic, 
P is the total power of the channels, 
d, is the noise power of the second order 
products within a band width of 
I ke/s, 
N is the number of channels and 
q is the power of the test signal referred 


to | milliwatt. 


picowatts 
non-linear phase-characteristic¢s 
in IF’-sections 


5 eee 
| | 
42 picowatts 
non-linear phase charac- 
teristics in RF-section 


The intermodulation noise power within 
one channel is proportional to the total 
harmonic power, to the density function of 
the second-order products, and to the ratio 
of the band widths. 

The amplitude-density function deter- 
mines the noise power within 1 e¢/s band- 
width. Harmonic distortion and linearity 
are interrelated, in case of quadratic charac- 


teristics, by the following equation : 


—m)* 
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where 


K, is the coefficient of the second order 
harmonic and 

m is the linearity defined by the for- 
mula 

—_ &min 


Sinadx 


me 


where gm is the mutual conductance. 

With quadratic characteristics assumed, 
there is a difference of about one order bhe- 
tween distortion coefficient and linearity. 

Noises due to non-linear phase charac- 
teristics are interrelated with the fluctuation 
of the group delay time by the following 
equations obtained by considerations similar 
to those applied above. 


r) = 1 
@ ae: 


In this equation the limit value of the 
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harmonic distortion is expressed as a func- 
tion of the drive and of the noise tolerated 
within one channel. Between the fluctuation 
of the time-delay and the limit value of the 


harmonic distortion the following equaiion 

holds : 

To 2 
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where 


t, is the fluctuation of the delay time, 
B is the IF band width, 
w isthe frequency of the channel 
under consideration and 
Aw is the frequency deviation. 


Since regarding noise problems the highest 
channel is operated under the most unfavour- 
able conditions, @ is the highest channel 
frequency. 

The reflection in the feed line and the 
noise power caused by the same are related 
in case of a long feed line by the following 
equation : 


NE i 2 ‘ 
d, = ae (az W (1AM) 


where 


u is the amplitude of the reflected wave, 

f is the frequency of the channel under 
consideration, 

t is the delay time of the one-way 
feed line, 

Ay is the peak frequency deviation and 

® is the echo-distortion function, 


The dependance on delay time t and 
peak deviation Ay can be given graphically. 

Using the above relations noise powers 
can be computed from the results obtained 
by measuring the non-linear properties of 
the amplitude and phase characteristics, and 
can be compared to the admissible values. It 
should be noted, however, that results ob- 
tained by linearity measurement calculations 
based on the same are only of informatory 
character. Intermodulation noises are deter- 
The 
telephone channels are replaced by a con- 
tinuous noise spectrum of appropriate band 
width. The band corresponding to the chan- 
nel under test is cut out from the noise 


mined by a special measurements : 


spectrum, and the noise power measured 
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within this band is characteristic to the inter- 
modulation distortion of the system. 

The results of the measurements 
have great influence to the rearrangement of 
the planned distribution of the admissible 
noise power. The final distribution will greatly 
depend on realization possibilities of the 
different components. Some types of modul- 
ator circuits, for instance, may cause distor- 
tions far below the admissible value. The 
noise power saved in this way can be used 


may 


to reduce requirements to be put to compo- 
nents in which the compensation of noise 
is involved with more difficulties. 

The C. C.1. T. T. recommendations con- 
cerning the hypothetical reference circuit 
constitute very high requirements regarding 
different sources of noise. 

Results of computations using the above 
relations and made for a 600-channel system 
show that in a band width of +10 Mes 
the admissible value of the delay time fluc- 
tuation is 1 nanosecond and that of the ampli- 
tude fluctuation 0,1 decibel: the reflection 
caused by the feed line should be less than 
1,05. The most considerable one of the RF- 
noises is the distant-field interference arising 
from the back radiation of the antenna. 
From this point of view, the optimum value 
corresponding to about 65 decibles can be 
obtained by the hog-horn type radiator. 
Noises due to this radiation can be compens- 
ated only by improving the signal to therm- 
al noise ratio, for the low value of the ad- 
missible intermodulation noise results in ex- 
tremely high requirements concerning each 
component part of the system. Jf the number 
of the channels to be transmitted by the 
system is raised above 600, the requirements 
outlined above are to be satisfied with a 
broader band width and with a higher drive. 
Owing to this fact the requirements concern- 
ing the individual components of the system 
are still higher. 


Miaiszaki felelés: Farkas Sandor 


A kézirat nyomdaba érkezett: 1958. IX. 16. — Terjedelem: 10,50 (A/5) iv, 56 abra, 1 melléklet 


Akadémiai Nyomda Budapest, V. Gerléezy u, — 46831/58 — Felelés vezeté: Bernat Gyirgy 


- CB and LB Sw 


ttchboards . Auto- 


: Microwave 
wave and Medium 
Mobile and P ; . Transmission 
Measuring 


: High Frequency Gene- 
rators 


BUDAVOX 


Budapest Telecommunication Company 
Budapest, y TT, Tandcs krt, 3Ja 
Telephone : 426-549 
Letters ; 

Budapest 62 p. O. B. 267 
Cables : 

Budavox Budapest 
Pelex : 672 
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RAILWAY CARRIAGE MANUFACTURERS 
and MECHANICAL ENGINEERS 
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Railcars and railcar-trains with diesel- 
mechanical, diesel-electric and diesel- 
hydraulic power transmission 
Passenger cars 
Special purpose railway vehicles 
Railway Diesel engines syst. Ganz ga rassik 
with output from 120 to 2000 HP 


Water turbines, Francis, Kaplan, Pelton type 
Hydro-electric power stations 
Micro-turbines, Ganz Mignon and Banki type 
Water pumps and pumping stations 


Letters: BUDAPEST, 10. P.O.B. 47 
Cables: GANZCOM, Budapest Telex: 532 
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oA Budapesti Muszaki Egyetem 
_. Periodica Polytechnica cimen ide- 
gen nyelvi tudomanyos folydiratot inditott. 


A folyoirat harom sorozatban — vegye- 
~ szeti, villamossagi, valamint gépész- es _ 
-altalanos mérnéki sorozatban — jelenik — 


-meg, évente negyszer, sorozatonkent egy- 
egy kétetben. Az egyes kdtetek terjedelme 
14—18 iv. 


A Periodica Polytechnicéban meg- 


- jelend tanulmanyok szerz0i az Egyetem 
tanari karabol és tudomanyos dolgozoibdl 
keriiinek ki. Fészerkeszt6 Dr. Cstros 
Zoltan egyetemi tanar, akadémikus. 

A folydirat eldfizetesi ara sorozatonkent 
és kdtetenként 60.— Ft. Megrendelheto 
az Akadémiai Kiadonal (Budapest 62, 
Postafiok 440. NB. egyszamlaszam: 09- 
915-111-44), a kiilfold szamara pedig a 
Kultura Konyv és Hirlap Kiilkereske- 
delmi_Vallalatnal (Bp. 62, Postafiok 149. 
NB. egyszamlaszam : 43-790-057-181): 
illetve a vallalat kulfoldi képviseleteinel es 
bizomanyosainal. 
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